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FOREWORD 

The research described herein was conducted by Willow Run Laboratories, a 
unit of The University of Michigan's Institute of Science and Technology, under 
Contract No. NAS 3-13037, sponsored by Lewis Research Center, National Aero- 
nautics and Space Administration. Project Manager for this contract was D. R. 
Buchele. F. S. Simmons was Principal Investigator. This final report covers the 
period from June 1969 through April 1971. The Willow Run Laboratories' report 
number is 30410-19-F. 
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ABSTRACT 

An investigation was conducted to explore the use of infrared and ultraviolet 
emission -absorption spectroscopy for determination of temperature profiles in 
flames. Spectral radiances and absorptances were measured in the 2.7-pm 
band and the 3064 -A OH band in H 2 /O 2 flames for several temperature profiles 
which were directly measured by a sodium line -reversal technique. The tempera- 
ture profiles, determined by inversion of the infrared and ultraviolet spectra, 
showed an average disagreement with line -reversal measurements of 50°K for the 
infrared and 200°K for the ultraviolet at a temperature of 2600°K. The reasons 
for these discrepancies are discussed in some detail. 
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MEASUREMENT OF TEMPERATURE PROFILES IN FLAMES 
BY EMISSION-ABSORPTION SPECTROSCOPY 

SUMMARY 


An experimental study was conducted to explore the use of infrared and ultra- 
violet emission -absorption spectroscopy in the 2. 7 -pm and 3064 -A bands of H2O 
and OH respectively, for the determination of temperature profiles in nonuniform 
flames . 

A rectangular multiple -diffusion burner with a number of unique features was 
designed and constructed for this study. A single manifold supplied oxidizer gas to 
the entire flame; the fuel manifold was divided into seven, isolated compartments. 
This permitted the establishment of temperature profiles along the length of the 
burner by throttling the fuel flow to each burner segment; thus, the local fuel/ 
oxidizer ratio was varied. To color only the central portion of flame, a method 
was developed to introduce sodium into the fuel gas supplying the central portion of 
flame, thus permitting a direct determination of the temperature profile by line- 
reversal measurements . 

A variation of the traditional line -reversal method was introduced for the 
temperature -profile measurement. The brightness temperature of a carbon arc 
was varied by an optical attenuator (designed and constructed for this purpose) 
consisting of a pair of counter-rotating, circular-disc, neutral density filters. 

With the attenuator set for the desired temperature, the line reversal was obtained 
by adjustment of the fuel -flow control valve, always keeping the mixture -ratio on 
the fuel-rich side of stoichiometric . 

The burner was successfully operated with H2/O2, CH4/O2. and CO/O2 as 
reactants . It was also operated in an isothermal mode, in which a central H2/O2 
flame was shielded by CO/O2 flames at either end to obviate the cooling and dilu- 
tion of the combustion product, H2O, mixing with ambient air. A 1.8-m Rowland 
circle spectrograph was converted into a scanning spectrometer by incorporation 
of an exit-slit assembly, which was designed and constructed for this study. A 
calibration procedure, which incorporated an NBS -certified tungsten lamp and a 
xenon arc, was developed for the reduction of the ultraviolet emission spectra into 
absolute spectral radiances. 

Spectral radiances and absorptances were measured in the infrared and ultra- 
violet for four temperature profiles in H 2 /0 2 flames: linear (viewed from either 
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direction); triangular, with a hot center; and triangular, with a cool center . The 
data are presented in complete detail, along with the sodium line -reversal tempera- 
tures . 


The molecular band model used for the inversion of the infrared spectra was 
modified to handle more accurately nonisothermal radiation calculations for hot, 
gaseous sources exhibiting large gradients in temperature. The improved band 
model was checked by a comparison of predicted and measured values for a small, 
methane -oxygen flame viewed through 10 m of ambient air. 

A computer code was written for use in inverting the OH spectra. It incorpo- 
rated a parametric approach similar to that used for the infrared. However, in- 
stead of using a band model, the spectral radiances were calculated for isolated 
spectral lines by integrating the exact equation of transfer over both the optical 
path and the contour of the spectral line. An analysis was made to determine the 
relative importance of pressure and Doppler broadening; the results indicated that 
the former dominates for the optical depths encountered in this study. 

The parametric method of inverting infrared emission and absorption spectra 
has been demonstrated in application to flames . Although the temperature profiles 
so determined showed disagreement with the profiles determined from the line- 
reversal observations, the disagreement is only nominal and the approach appears 
to be basically valid. In particular, the results indicate that a wide range of temp- 
eratures can be accommodated in such measurements . On the other hand, the 
parametric method of inverting ultraviolet emission and absorption spectra appeared, 
at least in this experiment, to be limited to providing information only over a narrow 
range of temperatures because of the steep variation of OH concentration with temp- 
eratures, especially in the region close to stoichiometric. Furthermore, the OH 
temperatures obtained were somewhat higher than the infrared and sodium line temp- 
eratures . This is tentatively attributed to OH chemiluminescence in and near the 
reaction zone. Further studies are required to delineate the regions and conditions 
in which OH chemiluminescence might be expected to affect emission absorption 
temperatures . 


2 



WILLOW RUN LABORATORIES 


1 

INTRODUCTION 

In a previous investigation at this laboratory [l], infrared emission -absorption 
spectroscopy was studied as a means for determining temperature profiles in inhomo- 
geneous hot gases, such as the exhaust jets of advanced air-breathing engines. A 
computer code, based on a generalized form of the statistical band model, was written 
for the calculation of molecular radiative transfer in diatomic and polyatomic gases. 
This code was applied to the determination of temperature profiles in hot gases by a 
procedure based on prior knowledge of the profile shape. This procedure involved 
the fitting of band -model predictions to observed spectra. For that previous study, - 
spectral radiances and absorptances were measured in the 2- to 3 -pm region for 
various pressures and temperature profiles along a 60 -cm path in furnace -heated 
samples of, H 2 O, CO 2 , HF, and H 2 O, CO 2 , N 2 mixtures . In all cases, the tempera- 
ture profiles. so obtained differed by an average value of 25°K from direct measure- 
ments by thermocouples on the sample cell . r 

In .the present study, a continuation of this research, emission -absorption ~ 
spectroscopy is applied to the determination of a higher temperature profile in an 
inhomogeneous flame to provide a more realistic simulation of actual applications . 

For this purpose, both infrared measurements in the 2.7-pm H 2 O band and ultra- 
violet measurements in the 3064 -A OH band were investigated. 
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2 

RADIATIVE -TRANSFER FORMULATIONS 


2.1. BAND -MODEL 

The band-model formulation used in this study is based on Eq. (5) of Ref. [ 2 ]; 
the rationale in the development of this model is reviewed in Ref. [l] . Expressed in 
terms of the average transmittance (f) in a spectral bandpass encompassing a number 
of lines, the average spectral radiance at the observer from a given optical path (E^ ) 
is given in approximate form by 

n 

A = 1 L ;< T i )<Af) i <« 

i=l 

in which L* is the Planck function of the temperature, T i . For this approximation, 
the optical path has been divided into n layers . The average spectral absorptance 
(a) of the optical path is given by 

n 

Oi(u) = 1 - ^ (AT). (2) 

i=l 

where v is the spectral frequency. The apparent (the indicated emission -absorption) 
temperature (T*) is given by 


L ^ T *> 


(3) 


The spectral transmittance difference, (AT) i = Tj-j — Tj, for a single gas is obtained 
from T j = exp [ -/? e f(Xj)], where p is the line -overlap parameter , and for a mixture 
of m absorbing gases from the relation 

m 


T 


n 



The function f(x) is the Ladenburg-Reiche expression or a simple algebraic approxi- 
mation of it [3] . The argument, x, is the dimensionless optical depth, defined by 

1 0<T.) - 


X. 

1 


= 0 


-1 


L k(T j>[ 

j=i 


T 


x ] 


(AX) 


J 


(4) 
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where k(Tj) = the local value of the first band -model parameter, the average absorp- 
tion coefficient (which can be identified with the ratio of average line 
strength to spacing) 

£(Tj) = the local value of the second parameter (identifiable as the ratio of 

average line width to spacing and proportional to the effective broaden- 
ing pressure) 

/Be = die absorption -coefficient weighted average over the entire path 

(AX)j = the optical depth increment 

The band -model parameters were extracted from the General Dynamics tabulations [4] 
and de’normalized as indicated in [1] .* The optical depth increment is given by (AX)j = 
Pj(AL)j, where Pj is the local value of the partial pressure of the radiating species and 
(AL)j is the geometrical path increment. The exponent V provides for the interpolation 
between the nearly -weak and nearly -strong line approximations . For calculations which 
involve only hot gases and in which the temperature gradient is not too severe, the em- 
pirically derived specifications x^ < 1, V = 0; 1 £ ^ 10, rj = (x^ - 1 )/9 ; x^ > 10, 

V = 1, where x^ is the optical depth for the entire path, are satisfactory. However, for 
calculations involving a large variation in temperature, a more complex interpolation 
procedure is necessary. Such a procedure has been developed and is described in Appen 
dix I. In Appendix II, some comparisons of the band-model predictions with hot gas 
spectra are shown. A particular application, the calculation of net radiation from hot 
combustion products viewed through an absorbing atmospheric path, is described in 
Appendix III . These changes have improved the general performance of the band model 
[l] . Appendix IV includes a discussion of some adjustments of the band-model param- 
eters for low temperatures . 

2.2. ISOLATED LINES 

Because of the simpler band structure of diatomic molecules, the rotational lines 
of OH were well separated and easily resolved by the instrumentation available. There- 
fore, like the infrared data for HF [1], the ultraviolet emission -absorption data could be 
handled on a line-by-line basis . However, approximate expressions for integrated line 
radiances and equivalent widths were avoided. Instead, these quantities were calculated 
using the exact equations of radiative transfer for a nonscattering source in local thermo 
dynamic equilibrium: 



*In Appendix C, p. 92 of [l], the denormalization. factor should read [273/T], not 
y 273 /tT 
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OL(V) 


1 - exp 



k(y,X) dX 


where k 
X 


the absorption coefficient 
the optical depth coordinate 


(5a) 


The integrated radiance and equivalent width are defined by 

( 6 ) 


(7) 

respectively, and the apparent temperature by 

L*(T*) = | (8) 

where W = the equivalent width . Equations (5) through (8) were programmed for nu- 
merical computations on The University of Michigan's IBM 360 system. These equations 
require that transition probabilities or oscillator strengths be transformed into the 
form of the line strengths used in the infrared. 

0 o 

The 3064 -A OH band results from transitions between the ir~ (lower) and the 
2 (upper) electronic states . The transitions in this band near 3064 A correspond 
to the v = v =0 band, where single and double primes indicate upper and lower 
levels specified by the vibrational quantum number, v. The rotational structure of 
this band has been described by Dieke and Crosswhite [6] . The line strength (S) for 
a particular line within the 3064 -A band is given by 


L 


00 

= J L 

J v 
0 


dv 


and 


W = a( v ) d„ 
J 0 


S 



k(^) dv 

tJ 

0 


(9) 


where k is the spectral absorption coefficient of Eq. (5) in cm’^ atm - -*-, and V is the 
wave number in cm"^ . The dependence upon the path has been suppressed. J / and ]" 
indicate the total angular momentum quantum numbers of the upper and lower rota- 
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tional states . The line strength is related to the oscillator strength fj"j' by [7]: 

2 


7re 


S j'Y . - 


( 10 ) 


me 


q 

where Nj" is the population of the lower level in cm , m and e are the electronic 
mass and charge; and c is the velocity of light. Because of the small population of 
the upper state, contributions to the line strength from emission by molecules in the 
upper state are neglected in this definition. The oscillator strength and rotational 
matrix element are related by 




x 


J 


R // / 

J J 


( 11 ) 


J (2J ' + 1) 

where Rj" t' is the relative magnitude of the matrix element, F 7 is an empirically 

determined constant, and T t " t 7 is a vibration -rotation interaction factor. For the 

J J - _ _ . r- . - _ „ A . _/ 

was 


tabulations of Rj" t' of Dielce and Crosswhite [6], a value of 2.0 x 10 4 for F * 


determined by Golden et al . [7], and values of T t " ' were computed by Learner [8] 
A combination of Eqs . (10) and (11) yields J 


T // / R 


N, 


Yy 


- M .1 J .1 J 


/ / 


(2 f 77 + 1)N 


(^) 


me 


( 12 ) 


where N is the total number density of OH molecules. With the prime notation omit- 
ted for convenience, the population factor Nj/N is expressed by the product 


N. 


N 



-E / K T 
e 



-E AT. 
v ) 



-E /KT 
r N 


(13) 


where E e , E v , E r are the energies for the electronic, vibrational, and rotational 
states respectively; g e , g v , g r are the corresponding statistical weights; and Q e , 
Qy, Q r are the electronic, vibrational, and rotational partition functions. The term 
in Eq. (13) corresponding to the electronic state is essentially unity for all tempera- 
tures of interest here . 


The vibrational energy levels were computed from the relation 
E v (v) = co(v + 1/2) -x(v + 1/2) 2 


(14) 
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where CO = 3735.21 cm'-*- and x = 82.21 cm‘1 [6]. The partition function, Q y , was 
computed by a summing of the exponential factors, exp [-E y (v)/K T] over v until con- 
vergence was obtained. Less than 10 terms were needed for convergence at 3000°K. 


The rotational state of an OH molecule is not specified until both K, essentially 
the rotational quantum number when electron spin is neglected, and J, the total angular 
momentum (J = K ± 1/2), are specified. Since in this study, only the R .2 branch of the 
0-0 band of the to ^2 transition is considered, by the Dieke and Crosswhite conven- 
tion, the relative population of molecules in the "2” state must be computed; that is, 
those molecules having not only a given value of K but a value of J = K - 1/2 as well . 
(In Dieke and Crosswhite, the line is designated by its K value; hence the R^l line has 
J = 1/2, K = 1.) Reference [6] gives expressions for the rotational energies of the 
"1" and "2" state of the lower (^tt) electronic state as follows: 


E (K) = B {(K + l) 2 - 1 - \J 4(K + l) 2 + A(A - 4) } - DK 2 (K + 1)' 
"1” state (J = K + 1/2) 


(15) 


E r (K) = B {k 2 - 1 + y 4K 2 + A(A - 4) } - DK 2 (K + 1)' 


2" state (J = K - 1/2) 


These expressions were used to compute the rotational energies with the following 
values for the rotational constants (B, D, A) appearing in Equation (15) [6]; 

B = 18 . 515 cm**- 
D = 0.00187 cm' 1 
A = -7.547 


The statistical weight is 2J + 1 for either of the two rotational states which corre- 
sponds to a given value of J. The rotational partition function was computed by eval- 
uation of the following sum: 


K 


max 


E r (K) 


Qr = ^ |(2K + 2) exp - ^ ^ — J + 2K exp 


-E r (K) 
2 

k T 


(16) 


K=0 


where k is Boltzman's constant. The sum was continued until convergence occurred 
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(i.e., to about K max = 40). For K > 100, the expressions for the rotational energy 
Equation (15) began to diverge; however, convergence of the sum in Equation (16) was 
obtained much before this occurred. Line strengths computed from Equations (12-16) 
for several lines at various temperatures are shown in Table 1 . 


Line widths (-y) for OH were taken from the recent work of Nadler and Kaskan 
[9] . Their resu lts are presented in terms of the ratio of Lorentz to Doppler half- 
width, a = vLnz yl/YD> typical values are shown in Table 2 for a pressure of one 
atmosphere and a temperature of 1485°K. In all cases, a > 0.2 so that, in accord- 
ance with the analysis by Yamada [10], the curve of growth of a Voigt line is indis- 
tinguishable from that of a Lorentz line over the entire range of optical depths . Hence, 
Lorentz line half -widths can be computed from values of a in Table 1 and adjusted to 
the local temperature and pressure, in accordance with the billiard-ball approxima- 
tion for collision broadening: 


V L (P. T) 


P 

7 L,o P 
o 



(17) 


where P Q = 1 atm, and T 0 = 1485°K. Table 1 lists values of -y l, for a number of lines 
at various temperatures . The values of a determined by Nadler and Kaskan are ques - 
tionable . Their work shows a sharp increase in a for large J values which is difficult 
to justify physically. On another contract, Anderson's theory was used to make calcu- 
lations of the line width of another diatomic molecule, HF, which has a similar rota- 
tional constant, B. The results indicate that, if anything, the line widths would be 
expected to decrease with increasing J value and, at high temperatures, would tend to 
uniformity over J. Moreover, there is a large amount of scatter in the results of 
Nadler and Kaskan, and their broadening gas (N 2 + T^O) was not the same as that of 
the present experiment (H 2 O). These factors indicate that a uniform value of a be- 
tween 0.3 and 0.5 corresponding to 1485°K for all lines might be more realistic for 
these experiments than the values given in Table 2. 


It is possible to determine a rotational temperature of the lower state from the 
rotational structure of the absorption spectrum from an isothermal flame. In this pro- 
cedure, the equivalent width of a single isolated line is given by 


W 



VV 

2ny 


X 


(18) 


where W is the equivalent width of the line having strength Sj" y' , and half -width y, 
and f is again the Ladenburg-Reiche function. X is the physical optical depth, F(OH)x 1 
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TABLE 1. STRENGTHS AND LORENTZ HALF-WIDTHS* 
OF SEVERAL OH LINES 



_ 1 


Strength 
(cm"2 atm' 1 ) 



Half Width [10] 
(cm -1 atm' 1 ) 


Line 

Vq cm 

1600 

2000 

2400 

2800K 

1600 

2000 

2400 

2800K 

V 

32415.51 

90.9 

71.2 

57.3 

47.0 

0.060 

0.054 

0.049 

0.045 

V 

32489.49 

257.0 

207.0 

170.0 

141.0 

0.045 

0.040 

0.037 

0.034 

V 

32517.58 

324.0 

268.0 

223.0 

183.0 

0.040 

0.036 

0.033 

0.030 

V 

32572.59 

324.0 

317.0 

286.0 

254.0 

0.040 

0.036 

0.033 

0.030 

V 4 

32551.72 

55.0 

84.0 

105.3 

118.7 

0.057 

0.051 

0.047 

0.043 

H 2 15 

32531.76 

36.0 

60.4 

80.8 

95.3 

0.057 

0.051 

0.047 

0.043 

v° 

* 

32363.49 
OH broadened 

| 2.63 7.72 

by H 2 0 + N 2 

15.0 

23.0 

0.095 

0.085 

0.078 

0.072 


TABLE 2. RATIO OF LORENTZ AND 
DOPPLER HALF-WIDTHS FOR SEV- 
ERAL OH LINES AT ONE ATMO- 
SPHERE AND 1485°K [10] 


Line 

A 

V 

0.50 

R 2 3 • 

0.37 

V 

0.33 

V 

0.32 

R 2 14 

0.46 

R 2 15 

0.46 

R 2 20 

0.70 
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where P(OH) is the partial pressure of OH, and 1 is the path length in cm reduced to 
standard temperature and pressure . The line strength St' ' y is a function of the rota- 
tional quantum number of the lower state, ]" , as given by Equations (12-16). Sj'U' 
can be expressed in the following form, in which the J dependence of Sj"j' appears* 
explicitly, but the temperature variation of S is combined into one common factor, 

K(T): 



// / 


J J 


K(T)'Tj"j/Rj'/ / exp 



(19) 


After substitution of this expression into Eq. (18) and after the inverse Ladenburg- 
Reiche function and the natural logarithm of both sides are taken, the following ex- 
pression is obtained: 



.-1 


{ 


(») 

' 27ry / 


1 2ny 


TV V j' 


-j2n [ K(T) - X] 


-E^") 

7t 


( 20 ) 


In a given measurement, 4.n [K(T) -X] is a constant for all the rotational lines in a 
band. Thus, if the first term on the left is plotted versus ~E r (J )/k for a number 
of rotational lines, the result should be a straight line the slope of which is 1/T 
where T is a rotational temperature characterizing the distribution of rotational 
states in the lower electronic state. In general, if this technique is to be applied, 
values for the line half -widths are required, and these are not accurately known. 

If the lines, as observed, are optically thin (in the region of linear growth), then 
the function f '•*-(W/27rY)27ry reduces simply to W and the line half -width no longer 
appears . Unfortunately, in these experiments, the lines were never in the linear 
growth region. 


; 
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3 

EXPERIMENTAL ARRANGEMENTS 

3 . 1. BURNER CONFIGURATION 

The nonisothermal flames studied in this program were created with a seven - 
section multiple diffusion burner schematically illustrated in Fig. 1 . The operation 
of this burner was similar in principle to one used in previous studies of HF flames 
[11] . The fuel gas passed through a set of copper tubes from seven compartments 
in the manifold; the flow into each was controlled by a separate valve. The oxidizer 
gas flowed between the tubes from a single manifold. Combustion occurred above 
the tube -end surface in a multiplicity of small diffusion flames; the combustion pro- 
ducts rapidly mixed to approximate conditions of a premixed flame. The temperature 
profiles were created by a throttling of the fuel gas entering each compartment. This 
caused a variation in the fuel -oxidizer ratio of the combustion above that segment. 

The controls were set to produce a desired temperature profile, as indicated by the 
sodium line -reversal temperature measurement described below . The top and side 
surfaces of the burner were cooled by water flowing through tubes of rectangular 
cross-section, which were soldered in place. 

The burner was constructed so that sodium could be introduced into the fuel 
gas flowing through the central tube in each segment (see Fig. 2). A line from the 
burner compartment passed fuel gas through the sodium vaporizer, then back to the 
compartment into the central tube of that segment of the burner. The vaporizer con- 
sisted of an electrically heated graphite boat containing a charge of sodium iodide. 
When the current through the boat stabilized, sodium -iodide smoke could be produced 
at a constant rate and would be picked up by the flowing fuel gas . Thus, each seg- 
ment of the burner could be individually colored by sodium for the line -reversal 
measurements which were used to set the temperature profiles . 

3.2. INFRARED INSTRUMENTATION 

The seven-segment burner was installed in the sample beam of a double-beam, 
infrared spectrometer shown schematically in Fig. 3. The configuration and opera- 
tion of this spectrometer were described in detail previously [l] . However, a change 
in the intensity calibration procedure was introduced in the present study. A separate 
blackbody was placed in the sample beam so that the effective infrared radiance temp- 
erature of the blackbody built into the spectrometer could be directly determined. 
Thus, the differences in mirror reflectances between 0.65 pm, the effective wave- 
length of the optical pyrometer previously used for calibration, and 2.7 pm, the wave- 
length of the infrared measurements, were completely accounted for . In a series of 
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FIGURE 1. SEVEN-SECTION, RECTANGULAR MULTIPLE -DIFFUSION BURNER 
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FIGURE 2. SODIUM -INJECTION SCHEME 
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FIGURE 3. EXPERIMENTAL ARRANGEMENT 
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measurements j the optical pyrometer was used to measure the temperature of the 
external blackbody and to monitor the internal blackbody . 

3.3. ULTRAVIOLET INSTRUMENTATION 

The ultraviolet measurements were made with a 1.5-m Rowland circle spectro- 
graph constructed some years ago at The University of Michigan’s Physics Department. 
The grating (an original, ruled by Rowland about 40 years ago) had a ruled area 90 X 55 
mm in size with 400 grooves/mm . The photographic plate holder was replaced by a 
movable exit-slit assembly constructed for this study, and a 1P28A photomultiplier 
tube was used as the detector . The grating was blazed for wavelengths in the visible; 
accordingly, the measurements of the 3064- A OH band were made in the second order . 

For the emission measurements, a 90 Hz chopper was positioned in front of the 
entrance slit; a Princeton Applied Research Model HR -8 Lock-in amplifier was used to 
amplify and rectify the photomultiplier signal for display on the strip-chart recorder. 
The intensity calibration was made with a National Bureau of Standards radiance stand- 
ard lamp positioned on the far side of the flame. The tungsten ribbon filament was 
imaged in the center of the flame and reimaged onto the spectrometer entrance slit 
by a pair of 250-mm focal-length fused quartz lenses . A Corning #9863 filter was used 
to eliminate the visible radiation from the standard lamp; the elimination was complete 
to within the noise level of the resultant measurements, as verified by temporary place- 
ment of an ultraviolet absorbing filter in the beam. Thus, the calibration procedure 
required only one correction factor to account for the reflection losses of the quartz 
lens behind the flame . 

The measurements of spectral absorptance were made with a 150-W xenon-arc 
source, chopped at a frequency of 2040 Hz and positioned behind the flame . The higher 
chopping frequency was necessary to discriminate against fluctuations in the flame 
which generally exhibit an amplitude varying inversely with frequency [ 12] . 

3.4. SODIUM LINE -REVERSAL INSTRUMENTATION 

Sodium line -reversal temperatures were measured transverse to the longer 
dimension of the burner by the apparatus shown in Fig. 3. A Mole-Richardson Type 
2371 carbon arc was used for the reference source. The arc crater was imaged in 
the center of the flame and reimaged onto the slit of a divided -circle spectroscope 
by a pair of lenses . The apparent brightness temperature of the arc source could be 
continuously varied by an attenuator consisting of a pair of counter rotating, circular 
disc, neutral density filters calibrated by means of an optical pyrometer. 
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The temperature profile along the length of the burner was established by the 
segment -by-segment sodium line-reversal measurements made in the transverse 
direction. First, the attenuator was set to the desired brightness temperature; then 
the fuel -flow rate was adjusted (on the fuel -rich side of stoichiometric) until the line 
reversal was observed. 

Unfortunately, the rectangular cross -section cooling tubes on the top of the 
burner required the transverse line -reversal measurements to be made with the mean 
height of the optical beam approximately 1 cm above the burner -tube ends . Conse- 
quently, the radial diffusion of the sodium from the central tube reached the mixing 
region of the flame gases with the ambient air at about the height of the optical beam 
(see Fig. 4). To correct the resultant error caused by absorption in the cooler mix- 
ing region, a set of line-reversal measurements were made at various fuel -air ratios, 
viewing the sodium -colored segment of the burner in both the transverse and longitu- 
dinal directions . In the latter case, of course, the uncolored adjacent segments , set 
at the same temperature, shielded the central segment. Thus, a correction was ob- 
tained for the transverse line -reversal measurements with which the temperature 
profiles along the burner were set; the magnitude of this correction is shown in Fig. 5. 

The experiments proceeded in the following sequence . A temperature profile 
was created and measured by the line -reversal observations . A complete set of in- 
frared and ultraviolet emission, absorption, and calibration spectra were then obtained 
with the individual fuel -line valves remaining at their set positions . Between scans, 
the burner was shut off by the two main flow valves, which were always set in the 
full -open position during burning . Thus, there was no problem in reproducing identi- 
cal profiles for the emission and absorption scans within a set. 
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FIGURE 4. DIFFUSION OF SODIUM IN FLAME 
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FIGURE 5. CORRECTION FOR LINE -REVERSAL MEASUREMENTS 
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4 

EXPERIMENTAL DATA 


4.1. INFRARED SPECTRA 

As in the previous study, the infrared spectra were recorded both on a chart 
and on cards for automatic reduction; for the latter, the sampling rate during scan- 
ning was sufficient to provide excellent fidelity of the reduced spectra to the chart 
recording. Data were obtained in this way for five settings of the burner for hydrogen - 
oxygen combustion. Linear profiles were created with the high- (profile A) and low- 
profile B) temperature ends toward the monochromator. (Actually, a single profile 
was created, and the burner rotated 180°.) Symmetrical triangular profiles were set 
with high - (profile C) and low- (profile D) temperatures at the center segment. Fin- 
ally, for a system operation check, an isothermal flame was produced (profile E) by 
adjusting the interior segments of the burner to a constant temperature, using H 2 and 
0 2 > while the segments on each end were set to the same indicated temperature using 
CO and 02- Thus, the nonisothermal mixing region was eliminated for observations 
in the 2. 7 -pm region outside the CO 2 band, but within the H 2 O bands . The reduced 
spectra for the nonisothermal runs (radiance, absorptance, and apparent temperature) 
are presented in Figs. 6 through 9; the reduced spectra for the isothermal run with 
the shielded flame are presented in Fig. 10. The corresponding line -reversal temper- 
atures are shown in Figs. 11 through 14; the burner segment designated 1 is closest 
to the monochromator . 

In the isothermal case, the apparent temperature, T*, should correspond to 
an actual flame temperature and be independent of wavelength. Figure 10 shows this 
to be the case except in the region of the 2. 7 -pm C0 2 band where mixing-region effects 
are noted, and at the extreme wavelengths, where measurement errors become large. 
Ignoring these regions, the resultant average apparent temperature of 2630°K agrees 
with the corrected sodium line-reversal temperature, 2600 K, at which the flame was 
set (Fig. 15) to within the overall experimental uncertainty. The source of the spec- 
tral noise in the trace of apparent temperature for the isothermal case (Fig. 10) was 
not determined; presumably, it is the consequence of imperfect register between the 
emission and absorption spectra in the data reduction. 

The burner was also operated with methane and carbon monoxide as fuels . For 
CH4/O2 operation, data were recorded and reduced; however, the quality of the results 
was lessened by a subsequent discovery of a procedural error in the setting of the flow 
rates, and circumstances did not permit repetition of these runs. For CO/0 2 operation, 
the much lower flame speed caused the flames to lift off the burner at comparable flow 
rates . For the flame to seat on the burner, a very low flow rate was required, and the 
small flame height did not adequately accommodate the optical beam. No data were re- 
corded for this case. 
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(a) Spectral Radiance 



FREQUENCY (cm' 1 ) 


. (b) Spectral Absorptance 
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FIGURE 6. REDUCED DATA, PROFILE A: LINEAR (HOT THROUGH COLD) 
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(a) Spectral Radiance 



(b) Spectral Absorptance 



FREQUENCY (cm' 1 ) 


(c) Apparent Temperature 

FIGURE 7. REDUCED DATA, PROFILE B: LINEAR (COLD THROUGH HOT) 
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FREQUENCY (cm' 1 ) 


(a) Spectral Radiance 



FREQUENCY (cm' 1 ) 


(b) Spectral Absorptance 



. (c) Apparent Temperature 


FIGURE 8. REDUCED DATA, PROFILE C: TRIANGULAR (HOT CENTER) 
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FREQUENCY (cm" 1 ) 
b. Spectral Absorptance 



FIGURE 9. REDUCED DATA, PROFILE D: TRIANGULAR (COLD CENTER) 
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FIGURE 10. REDUCED DATA, PROFILE E: ISOTHERMAL 
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4.2. ULTRAVIOLET SPECTRA 

The instrumentation for the OH measurements could not easily record on 
punched cards . Therefore, the data were recorded only on a strip chart and the 
reduction performed by hand. However, since the lines were isolated and the slit 
function of the spectrograph approximately triangular, the integrated radiances and 
equivalent widths could be taken as proportional to the heights of the lines on the 
charts . Errors caused by this procedure would mostly cancel in the determination 
of apparent temperature (Eq. 8), provided that the apparent widths of the emission 
and absorption lines are the same; the remaining errors were of the order of the 
experimental uncertainty. A more detailed reduction of the absorption spectrum of 
the isothermal flame was performed for the rotational temperature analysis discuss- 
ed later. In this case, equivalent widths were obtained with a planimeter . 

y 

Figure 16 shows a portion of the absorption spectrum of an isothermal H 2 /O 2 
flame obtained under the best resolution possible with this instrument. The lines 
were identified by comparison with the spectra obtained by Kostkowski and Broida [13] . 
Unfortunately, the energy in emission was not sufficient to be measured with an accept- 
able signal -to -noise at this 10-pm slit opening. Accordingly, the sets of emission and 
absorption spectra were obtained with a 25 -pm slit opening. 

Emission and absorption in the R 2 branch of the v = 0 0 band of OH were 

measured for the five profiles in H 2 /O 2 flames described by Figs . 11 through 15, as 
they were for the infrared. The unreduced emission and absorption spectra, with the 
calibration traces superimposed, are presented in Appendix V. 

Reduced data, the total radiances, equivalent widths, and apparent temperatures 
obtained for seven well -isolated lines (R 2 L 3, 4, 7, 14, 15, and 20) are presented in 
Table 3 . 

4.3. ISOTHERMAL TEMPERATURES 

If the infrared and ultraviolet portions of the spectrum are to be useful together 
for the prediction of temperature profiles, the processes involved must be shown to 
yield identical temperatures in an isothermal flame. Such a requirement generally 
means that the processes must be in thermodynamic equilibrium. The emission -absorp 
tion temperatures measured for the isothermal flame in the infrared region agreed with 
the sodium line -reversal temperatures within the overall experimental accuracy. How- 
ever, as indicated in Table 3 and in Fig. 15, the emission -absorption temperatures 
measured for this profile in the ultraviolet are approximately 200°K higher, a discrep- 
ancy which cannot be explained by radiometric errors in the experimental procedure. 
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FIGURE 16. ULTRAVIOLET ABSORPTION SPECTRUM OF ISOTHERMAL FLAME 
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The uncertainty in the radiance calibration (about 25°K) is far too small to account 
for such a large discrepancy. Also, scattered light and stray light from the first 
order of the grating were known to be negligible . 

One is then tempted to ascribe the high emission levels to an excess population 
of ^2 + OH radicals . This would indicate the prescence of a departure from thermal 
equilibrium and, hence, chemiluminescence in the gases above the reaction zone. A 
number of authors [ 14-18] have observed exactly these effects to a very significant 
extent within and near the reaction zone of premixed flames at atmospheric pressure, 
presumably caused by radical recombination . Such luminescence might be a possible 
explanation for the high ultraviolet emission -absorption temperatures observed here. 

On the other hand, the reaction zone of the flame in the current experiment, 
in the center burner segment at least, appeared to be below the slit image . The 
reaction zone could be seen to begin about 1/16 in. above the burner and to extend 
as a rather diffuse layer to perhaps 5/32 in. above the end of the burner tubes . 

The image of the spectrograph and spectrometer slits was located, although im- 
precisely, above this reaction surface. We estimate that the lower boundary of the 
slit image was about 1/8 in . above the burner in all cases . The gas velocities were 
such that the gases took approximately 10 ~ 3 sec to travel from the bottom of the re- 
action zone to this lower boundary of the slit image. Most evidence indicates that at 
atmospheric pressure, the OH should have equilibrated by this time. The slit image 
was located over the center burner segment during all measurements and, because 
of the above reasoning, the gas measurements for this segment should have been made 
in a region where thermodynamic equilibrium was established. However, because of 
the divergence of the beam between the slit image and the objective lens 20 in . away 
from the slit image (Item 22 in Fig. 3), portions of the beam passed much closer than 
1/8 in . to the burner in the outer segments . We believe that this passage of a fraction 
of the beam very close to the burner surface is the best explanation for the high emission - 
absorption temperatures in the ultraviolet. When the beam passes within 1/8 in. of the 
burner, emission from the reaction zone, including OH chemiluminescene, will be 
sensed. Because of the imprecision in the slit image position, we estimate that as much 
as perhaps 10 percent of the beam volume was within or near the reaction zone. How- 
ever, considering the levels of shemiluminescence observed in [14-18], this percentage 
does not seem sufficient to explain the high emission values leading to the high emission - 
absorption temperatures . The emission values are such that about 3.5 times as much 
OH emission is observed as would be observed from the combustion products at equili- 
brium at a temperature of 2600°K. It is possible, that, considering the variety of re- 
sults observed in the entries of [ 14 - 18 ] and considering that the current situation is 
greatly different from any of the situations described in [14-18] , the region of chemi- 
luminescence might extend significantly above the apparent reaction zone so that much 
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more than 10% of the beam could observe OH chemiluminescence. 

The OH absorption values were affected much less by this beam divergence 
because, although the amount of OH in the reaction zone is above the equilibrium 
amount, it does not exceed the equilibrium amount by orders of magnitude, and most 
of it remains in the ground state. Conversely, the emission is governed by the popu- 
lation of the ^2 ^ state of OH, which is orders of magnitude greater in the reaction 
zone than its extremely small value at equilibrium at 2600°K. The infrared measure- 
ments are not severely affected by this beam divergence because the vibrational states 
of H 2 O remain near equilibrium even within the reaction zone, although the H 2 O con- 
centration is somewhat lower there . 

In other experiments where OH chemiluminescence was observed, the rota- 
tional temperature of both the ir and the 2 states of OH were near the flame temper- 
ature. To check to see that our measurements concur with the result, an analysis of 
the rotational structure of the isothermal flame (Profile E) spectra was performed in 
absorption with the use of Eq. (20). When the same line width for all rotational lines 
was used, temperatures in absorption near 2570°K were obtained for line widths of 
0.03 and 0.05 cm' 1 (corresponding to values of a equal to 0.3 and 0.5 at 1485°K). 

The results are shown in Fig. 17a and agree to within the experimental uncertainty 
with the infrared and sodium line -reversal temperatures. With the values for the line 
widths obtained from Nadler and Kaskan (shown in Table 2), an unusually low tempera- 
ture in absorption of 2260°K was obtained (Fig. 17b) which casts further doubt on the 
applicability of their results to this study. This analysis shows that, while a uniform 
change in the line half -width of all rotational lines appears to have little effect on the 
rotational temperature obtained, the assumed variation of line half -width with rotational 
quantum number appears to have a great effect . The different rotational temperatures 
obtained illustrate the difficulties in the performance of a rotational temperature analy- 
sis when the lines are in the region of square-root growth. 

Accordingly, we have ascribed the discrepancy between the ultraviolet emission - 
absorption temperatures and the infrared and sodium temperature to the inclusion of 
a region of OH chemiluminescence in part of the beam. It was originally hoped that 
we would be able to use both the infrared and ultraviolet measurements together in the 
inversion procedures, but the temperature discrepancies made such a procedure im- 
practical. In the inversion procedures, the ultraviolet and infrared results were each 
considered separately. 
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(a) With Constant Line Widths of 0.03 and 0.05 cm" 1 
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(b) With Line Widths Derived from Nadler and Kaskan [9] 


FIGURE 17. RESULTS OF ROTATIONAL TEMPERATURE ANALYSIS OF 
• ISOTHERMAL ABSORPTION SPECTRUM 
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5 

INVERSION PROCEDURES 

There are a number of distinct approaches that have variously been utilized 
for inversion of spectroscopic measurements of flames and hot gases. These could 
be categorized as being either formal or empirical. Each category encompasses a 
considerable variation in approach. To accomplish the inversion, a formal approach 
could be based on either a mathematical operation, such as the Abel transform [12] 
or on a numerical computation, in which a set of radiative -transfer equations are 
simultaneously solved. In the latter approach, a set of n equations (e.g. , Eq. 4), 
representing measurements of emission and absorption at n wavelength intervals, 
could, in principle, be solved to yield profiles in temperature and concentration char- 
acterized by a total of n parameters . The problems in performing such calculations 
and obtaining convergence with real data with finite uncertainties are discussed in the 
literature . 

Recently the empirical approach has been found to be more useful. This ap- 
proach utilizes auxiliary information to insure convergence in the inversion compu - 
tation. For example, in the inversion of spectral data obtained with satellite-borne 
sensors, profiles in atmospheric temperature and water -vapor concentration are 
rapidly obtained by use of the nominal profiles of a standard atmosphere as a starting 
point in the iteration [19]. Analogously, determination of temperature and concentra- 
tion profiles in gas dynamic systems such as combustion or propulsion devices is 
easier if we have prior knowledge of the general shape of the profiles to be expected . 

This information is usually available from theoretical gas dynamic considerations . 

Thus, in the prior study of temperature -profile determination in hot furnace -heated 
gas samples [ 1 ], the inversions were based on the assumption of known profiles, 
relatively simple in shape, which would usually be characterized by one to three 
parameters . 

In the present study, the same basic approach was taken. However, in the case 
of flames, both the concentration and temperature profiles were unknown, so that the 
inversions were more difficult. For this reason, the temperature profiles studied 
were limited to simple shapes, linear and symmetrical triangular, which could be 
characterized by two parameters, the maximum and minimum temperatures. Further- 
more, because of the nearly linear relationships between the theoretical adiabatic 
flame temperature and the partial pressures of H 2 O (shown in Fig. 18), the same 
shape could be taken for the profiles in this quantity also. The theoretical temperatures 
and compositions for adiabatic H 2 0 , CH 4 -O 2 and CO-O 2 flames were obtained from 
the NASA thermochemistry computer -code [20] ; the values are presented in Tables 4-6. 
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TABLE 4. THERMODYNAMIC PROPERTIES OF ADIABATIC H 2 /0 2 FLAME 


Oxidizer to Fuel Ratio 

2.0 

4.0 

6.0 

8.0 

10.0 

Equivalence Ratio 

3.97 

1.984 

1.323 

0.992 

0.794 

P (atm) 

1.0 

1.0 

1.0 

1.0 

1.0 

T (°K) 

2065 

2837 

3047 

3078 ’ 

3054 

Enthalpy (cal/gm) 

0.0 

0.0 

0.0 

0.0 

0.0 

Entropy (cal/gm°K) 

8.446 

6.080 

4.978 

4.339 

3.921 

Average Molecular Weight 

6.041 

9.649 

12.56 

14.90 

16.82 

(91nV/ 31nP) x 

- 1.001 

- 1.023 

- 1.063 

- 1.077 

- 1.068 

(91nV/31nT)p 

1.015 

1.466 

2.246 

2.518 

2.363 

Specific Heat at Constant Pressure (cal/gm°K) 

1.661 

3.100 

4.819 

4.771 

3.912 

Specific Heat Ratio 

1.255 

1.137 

1.115 

1.112 

1.112 

Sonic Velocity (m/sec) 

1889 

1667 

1500 

1382 

1295 

Mole Fractions 






H 

0.0021 

0.0624 

0.0906 

0.0765 

0.0570 

H 2 

0.7461 

0.4558 

0.2507 

0.1486 

0.0952 

h 2 o 

0.2516 

0.4610 

0.5585 

0.5811 

0.5745 

o 

- 

0.0016 

0.0156 

0.0336 

0.0443 

OH 

0.0001 

0.0185 

0.0705 

0.1070 

0.1209 

0 0 

- 

0.0007 

0.0140 

0.0531 

0.1081 


36 



WILLOW RUN LABORATORIES 


TABLE 5. THERMODUNAMIC PROPERTIES OF ADIABATIC CH 4 /0 2 FLAME 


Oxidizer to Fuel Ratio 

1.0 

2.0 

3.0 

4.0 

5.0 

Equivalence Ratio 

3.99 

1.995 

1.330 

0.997 

0.798 

P (atm) 

1.0 

1.0 

1.0 

■ 1.0 

1.0 

T (°K) 

1039 

2618 

3031 

3053 

3018 

Enthalpy (cal/gm) 

-557.7 

-371.8 

-278.9 

-223.1 

-185.9 

Entropy (cal/gm°K) 

4.167 

3.883 

3.473 

3.193 

2.997 

Average Molecular Weight 

11.60 

15.80 

19.21 

21.50 

23.23 

OlnV/01nP) T 

-1.064 

-1.008 

-1.057 

-1.070 

-1.062 

0lnV/31nT)p 

1.954 

1.166 

2.148 

2.415 

: 2.282 

Specific Heat at Constant Pressure (cal/gm°K) 

3.262 

1.104 

2.992 

3.194 

2.773 

Specific Heat Ratio 

1.158 

1.173 

1.115 

1.110 

1.110 

Sonic Velocity (m/sec) 

928.7 

1271 

. 1209 

1145 

1095 

Mole Fractions 






CO 

0.2729 

0.2913 

0.2153 

0.1547 

0.1124 

co 2 

0.0204 

0.0371 

0.0840 

0.1134 

0.1289 

H 

- 

0.0245 

0.0645 

0.0499 

0.0351 

% 

0.6128 

0.3542 

0.1396 

0.0731 

0.0447 

h 2 o 

0.0371 

0.2880 

0.3953 

0.3913 

0.3708 

0 

- 

0.0002 

0.0179 

0.0392 

0.0479 

OH 

- 

0.0047 

0.0630 

0.0938 

0.0994 

o 2 

- 

0.0001 

0.0205 

0.0848 

0.1608 

c(s) 

0.0316 

- 




CH 4 

0.0251 

- 
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TABLE 6. THERMODYNAMIC PROPERTIES OF ADIABATIC CO/Og FLAME 


Oxidizer to Fuel Ratio 

0.150 

0.300 

2 

0.450 

0.572 

0.750 

Equivalence Ratio 

1.584 

1.311 

1.119 

0.999 

0.865 

P (atm) 

1.0 

1.0 

1.0 

1.0 

1.0 

T (°K) 

2198 

2915 

2975 

. 2977 

2958 

Enthalpy (cal/gm) 

-820.1 

-725.4 

-650.4 

-599.9 

-538.9 

Entropy (cal/ gm°K) 

2.083 

2.099 

2.094 

2.087 

2.076 

Average Molecular Weight 

32.21 

34.61 

35.16 

35.42, 

35.68 

(31nV/31nP) T 

- 1.000 

-1.032 

-1.052 

-1.055 

-1.053 

(31nV/31nT) p 

1.000 

1.722 

2.133 

2.204 

2.165 

Specific Heat at Constant Pressure (cal/gm°K) 

0.3195 

1.251 

1.722 

1.795 

1.744 

Specific Heat Ratio 

1.240 

1.116 

1.108 

1.107 

1.107 

Sonic Velocity (m/sec) 

838.6 

883.8 

882.8 

879.4 

873.4 

Mole Fractions 

CO 

0.7374 

0.5350 

0.4194 

0.3491 

0.2703 

co 2 

0.2626 

0.4155 

0.4463 

0.4553 

0.4576 

O 

- 

0.0154 

0.0329 

0.0409 

0.0462 

°2 

- 

0.0342 

0.1015 

0.1547 

0.2260 
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Since the H 2 -O 2 flames were found to be non -adiabatic, we have included Table 7, also 
obtained from Ref. 20, which gives the mole fraction of OH for such flames . 

Since there were four unknowns, at least four spectral data points were required 
to accomplish the inversion of the infrared spectral data. One choice would have been 
four values of the combined emission -absorption spectrum, expressed in the apparent 
temperature, T* . Such a choice would have preserved consistency with the previous 
study in which temperature profiles were deduced by matching curves of observed val- 
ues and band -model predictions of T* as a function of absorptance. In the present study, 
it was more practical to utilize the emission and absorption data separately at two wave- 
lengths representing fairly high and fairly low optical depths; the frequencies 3100 and 
3400 cm’ 1 were selected for this purpose. The values of the maximum and minimum 
temperature and water partial pressures were adjusted iteratively in calculations using 
Eqs . (1) and (3) until a best fit was obtained at the selected wavelengths in the emission 
and absorption spectra. (Actually, values of apparent temperature were used in place 
of radiance for convenience . ) The OH inversion proceeded analogously . Since the lines 
were well separated, Eqs . (6) and (7) were used for the theoretical predictions with the 
line strengths given in Table 1 . The line widths of Nadler and Kaskan (Table 1) were 
also used, except where noted . 

Initially, it was thought that the flames would be adiabatic, so that there would 
be a one-to-one correspondence between the local temperature and the mole fraction of 
OH according to Fig. 18 . However, the temperatures determined by the various methods 
for the isothermal profile, which was set close to a stoichiometric mixture ratio, are 
considerably lower than the theoretical values, indicating a departure from adiabaticity . 
This failure to achieve the theoretical flame temperatures indicated that it would no long- 
er be valid to use Fig. 18 to relate the temperature to the OH mole fractions, and hence 
that the OH -inversion would require four parameters also. 

In this case, however, a linear variation in OH pressure over the entire length 
of the burner would not be reasonable, as the degree of dissociation of H 2 O is a steep 
function of temperature . An alternative four -parameter iteration procedure would have 
involved compound linear variations as specified by the quantities Tmax- F(OH)max> 
T m in» and Xj^, where X m i n represents the distance along the burner at which the 
pressure of OH is taken to be zero, and T m in is the corresponding temperature . How- 
ever, a continuous variation in Xj^^ could not be accommodated because the computation 
was set up in terms of seven isothermal regions corresponding to the burner segments. 
Therefore, a different approach was taken to treat the OH spectra. 

To determine the OH concentration as a function of temperature, the absorption 
spectrum of the isothermal flame was examined in detail. The sodium line-reversal 
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TABLE 7. MOLE FRACTION OF OH FOR NON -ADIABATIC 
CONDITIONS, H 2 -0 2 COMBUSTION AT 1 ATM. 


O/F 

i 2.0 

4.0 

6.0 

8.0 

10.0 

T,°K 

• 





3100 

.01990 

.04600 

.079 20 

.10980 

.12760 

3000 

.01420 

.03370 

.06180 

.09250 

.11030 

2900 

.00980 

.02360 

.04560 

.07520 

.09260 

2800 

.00650 

.01580 

.03170 

.05910 

.07560 

2700 

.00410 

.01010 

.02090 

.04490 

.06010 

2600 

.00250 

.00620 

.01310 

.03290 

.04650 

2500 

.00150 

.00370 

.00780 

.02340 

.03500 

2000 

.00005 

.00011 

.00024 

.00250 

.00510 

1— * 
cn 
O 

o 

.00000 

.00000 

.00000 

.00009 

.00019 

T( adiabatic) 

2065 

2837 

3047 

3078 

3054 

OH (Mole Fraction) .0001 

.0185 

.0705 

.1070 

.1209 
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temperature, the infrared emission -absorption temperature, and the rotational tempera- 
ture obtained from an analysis of the relative line intensities all gave temperatures 
near 2600°K. Thus, it seemed reasonable to assume 2600°K for the vibration -rotation 
temperature of the flame. By using Eq. (18) with the line intensities of Table 1 and a 
temperature of 2600°K, an estimate was obtained of the OH concentration required to 
give the measured equivalent widths for the 7 lines considered. One estimate was ob- 
tained when the line widths of Nadler and Kaskan (Table 1) were used, and one estimate 
was obtained with a uniform line width of 0.04 cm'^. A value for the mole' fraction of 
OH of 0.028 was obtained in both cases. This value is considerably below the OH con- 
centration for a stoichiometric adiabatic flame and considerably above that for an adia- 
batic flame at 2600°K, further indication of the departure from adiabatic conditions . 
Reliable experimental oxidizer -fuel ratios could not be obtained, because the fuel con- 
trol valves could not be calibrated. However, since we were striving for near -maximum 
temperatures during the isothermal runs, the oxidizer -fuel ratio for these runs would 
have been slightly less than 8; the value corresponding to the maximum adiabatic temp- 
erature. The value obtained from the line strengths for the OH mole fraction, 0.028, 
compares favorably with the value of 0.033 in Table 7 for 2600°K and an oxidizer -fuel 
ratio of 8 . This agreement also indicates that the majority of the beam passes through 
a portion of the flame where OH concentration is near its equilibrium value . Table 8 
shows the observed equivalent widths for these lines and computed equivalent widths 
for two OH concentrations and for both choices of line width . 

The inversion procedure was applied in the ultraviolet to emission -absorption 
temperatures which are higher than the vibration -rotation temperature. Because temp- 
eratures could not be correlated with mixture ratios it was not possible to use Table 7 
to relate temperature to OH concentration, hence the adiabatic curve for OH concentra- 
tion versus temperature shown in Fig. 18 was modified as follows: To ensure that the 
OH concentration and emission -absorption temperature predicted by the inversion pro- 
cedure were consistent with the observed values in the isothermal flame, the theoretical 
curve of P(OH) versus T shown in Fig. 18 was shifted to the left to coincidewith the values 
T = 283 5°K and B(OH) = 0.028atm; i.e., the effect of flame cooling is assumed to be 
independent of the flame temperature and the discrepancy in OH temperature. This 
procedure yielded a new curve representing a one-to-one correspondence of OH partial 
pressure with ultraviolet emission -absorption temperature (Fig. 19). (In retrospect, 
a preferable procedure would have involved a set of high -resolution absorption spectra 
for isothermal shielded flames at various temperatures.) 

This inversion, then, proceeded like the inversion for the infrared data. A pair 
of lines R 2 1 and R 24 , were selected to represent large and small optical depths, re- 
spectively. The values of maximum and minimum temperatures (and the corresponding 
OH concentrations from Fig. 19) were iteratively adjusted in calculations in which Eqs.(6) 
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TABLE 8. EXPERIMENTAL AND THEORETICAL VALUES OF THE 
EQUIVALENT WIDTH OF OH LINES IN ISOTHERMAL 
FLAME AT 2600°K 



Observed 

Equivalent 

Width 


Calculated Equivalent Width (cm *) 


Line Half -Width = 0.04 cm * 

Line Half -Widths from Table 1 

Line 

(cm"l) 

P(OH) = 0.055 

P(OH) = 0.02 

P(OH) = 0.055 

P(OH) = 0.02 

V 

0.675 

0.887 

0.529 

0.965 

0.580 

R 2 3 

1.05 

1.55 

0.933 

1.46 

0.880 

r 2 4 

1.33 

1.77 

1.07 

1.57 

0.945 

R 2 7 

(1.46)* 

2 .04 

1.23 

1.82 

1.09 

R 2 !4 

0.863 

1.31 

0.785 

1.39 

0.843 

R 2 15 

0.797 

1.16 

0.695 

1.24 

0.733 

R 2 20 

0.374 

0.532 

0.309 

0.715 

0.400 


* Questionable value caused by slight overlap by adjacent line. 



TEMPERATURE (°K) 


FIGURE 19. MOLE FRACTION OF OH AS A 
FUNCTION OF TEMPERATURE 
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and (7) were used. Essentially two temperatures have appeared in the analysis to this 
point, the rotational -vibrational value measured by three methods and the electronic 
value measured by ultraviolet emission -absorption spectroscopy. In performing this 
iteration, one must now choose between the two values. If nonequilibrium effects are 
occurring, Eq. (5) is no longer valid, and the interpretation of the emission spectrum 
becomes more difficult. However, since we have ascribed the high ultraviolet temper- 
ature to the inclusion of only a small region where non -equilibrium effects are present, 
we have performed the inversion, using Eq. (5), based on the ultraviolet emission- 
absorption temperatures . This provides inversion results which agree with the meas - 
ured integrated radiances . 

By this means, we expected to obtain temperature profiles from a simple two- 
parameter fit of the integrated radiances for these two lines, letting the equivalent 
widths fall where they may. However, after a few iterations, it became obvious that 
because of the very steep gradient in OH partial pressure, no information could be 
obtained about the lower temperature regions in the flame . Test cases for each of the 
profiles, A through D, showed that a change of 100°K in the minimum temperature 
would yield a negligible change in the values of the integrated radiances and equivalent 
widths of the two lines. The results of these calculations, in effect, yielded only values 
for the maximum temperatures. Clearly, the experimental study of temperature meas- 
urements by OH emission and absorption indicate that they can contribute only to the 
determination of the maximum temperature. 

At the beginning of the present study, the intention was to utilize the results of 
the above four- and two -parameter fitting procedures as the starting point for a second 
approximation to the temperature profile. This approximation would have been charac- 
terized by 14 parameters, in correspondence with the 7 segments of the burner. For 
this purpose, 7 wavelengths in the 3100- to 3700-cm‘l region would have been selected 
for the infrared inversion, and 7 well isolated OH lines for the ultraviolet, as specified 
in Table 2. An iterative calculation, similar to that proposed by Smith [19], would 
have been employed. However, limitations in time and funds did not permit the inclu- 
sion of this analysis . 
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6 

TEMPERATURE -PROFILE DETERMINATIONS 


6.1. INFRARED 

The temperature profiles obtained from the inversions of the infrared spectra 
when the procedures detailed in the preceding section were used are shown in Figs . 11 
through 14. The results are in qualitative agreement with the temperature profiles ■ 
measured directly by the line -reversal apparatus. As the iterations proceeded, the 
convergence upon the temperature and concentration profiles representing the best fits 
to the observed spectra became quite slow. For expedience, the iterations were ter- 
minated when the calculated apparent temperatures and absorptances reasonably 
matched the observed values . The resultant best -fit conditions are summarized in 
Table 9. For profiles A through C, the observed and calculated values of T* and a 
agreed to within 25°K and 0.02 at the two selected frequencies . However, the conver- 
gence for profile D was much slower, and the calculations were terminated sooner. It 
was perfectly clear, however, that the results would not have been appreciably altered 
by further iterations . Therefore, the discrepancies in the temperature profiles as 
determined from the infrared spectra and the D-line observations remain to be explained. 
For the linear profiles, the greater differences occur at the lower temperature. (Recall 
that profiles A and B were achieved by merely rotating the burner 180° .) On the other 
hand, for the triangular profiles, C and D, the larger discrepancies occur at the higher 
temperatures . 

The reason for this effect is not at all obvious . However, there are several 
possible sources of error: 

(1) As mentioned in a preceding section, the D-line reversal measurement was 
affected by the addition of the cooling tubes to the top of the burner and the consequent 
raising of the optical path to a height at which some of the mixing region might have been 
colored with sodium. It was observed that the reversal did not occur consistently with 
perfect blending of the D -lines into the background continuum, indicating some variation 
in temperature along the path . 

(2) In the iterative spectral radiance and absorptance calculations made with 
the bandmodel, the mixing region at each end of the burner was not considered to 
avoid complicating the inversion procedure. Absorption in the cool mixing region 
would most affect calculations in which the higher temperatures are closer to the 
observer . It is felt that this probably accounts for the poor fit of the observed and 
calculated values of T* and OL shown in Table 9 for profile D . 
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(3) Another shortcoming in the calculation procedure involved mixing of flames 
from the adjacent burners . The calculations were made under the assumption that the 
optical path consisted of seven isothermal sections, each the width of a burner segment; 
inclusion of mixing effects between sections would have considerably complicated the 
analysis. During the iterative calculations, it was obvious that the absorptance at both 
frequencies, 3100 and 3400 cm’^, depended strongly on the average H 2 O concentration 
and only showed weak response to moderate variations in the concentration profile. This 
is consistent with the rather slow variation of H 2 O concentration with temperature, which, 
in turn, governs the variation in spectral absorption coefficient. However, the resultant 
H 2 O -concentration profiles shown in Table 10, although they are accurate to no more than 
10%, are clearly higher than the corresponding theoretical values shown in Fig. 18. This 
tends to substantiate the conclusion that the flames were not adiabatic. 

6.2. ULTRAVIOLET 

The maximum temperatures resulting from the inversions of the OH emission- 
absorption spectra are also shown on Figs. 11 through 14. Clearly, to display the com- 
plete profile of the final iteration would make no- sense as the uncertainty at the lower 
end would be enormous . However, the values used, as given in Table 10, are probably 
realistic, since the calculated and observed equivalent widths are in fair agreement. 

The values of the maximum temperatures so determined are considerably higher 
than those deduced from the infrared inversions, to a degree comparable to the differ- 
ences noted in the direct emission -absorption measurements (Fig. 15). This is by no 
means coincidental, since the OH pressure -temperature relationship is based on this 
very data . 
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TABLE 10. COMPARISON OF ULTRAVIOLET SPECTRAL DATA 
AND BEST-FIT CALCULATIONS 


P L W 

max 0 



T 

max 

T . * 
min 

max 
for OH 


(jLtW/cm^ 

-sr) 

(cm 

_1 ) 

Profile 

(°K) 

(°K) 

(atm) 

Line 

Observed Calculated 

Observed 

Calculated 

A 

2900 

1600 

0.043 

R„1 

0.82 

0.93 

0.306 

0.450 





2 

R 2 4 

1.50 

1.40 

0.678 

0.763 

B 

2860 

1600 

0.033 

V 

1.17 

1.07 

0.306 

0.409 





R 2 4 

1.70 

1.82 

0.678 

0.694 

C 

2940 

1700 

0.060 

V 

1.64 

1.70 

0.382 

0.488 





R 2 4 

2.73 

2.56 

0.736 

0.831 

D 

2940 

1700 

0.060 

V 

3.17 

3.03 

0.538 

0.648 





V 

4.60 

5.01 , 

0.983 

1.092 


* Highly uncertain 
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7 

CONCLUSIONS AND RECOMMENDATIONS 

As a result of the current experimental program, it appears that, in the infrared 
region at least, temperature profiles in flames can be measured by emission -absorption 
spectroscopy in the same manner as the hot gases of the previous study [Ref. l] . The 
uncertainties in the temperatures deduced during the current series of measurements 
were approximately 30° to 50°K while the temperature uncertainty in the earlier program 
was about 10°K. 

An error analysis was undertaken which considered the temperature differences 
(2600°K versus 1200°K) and absorption differences (~ 1 .0 versus ~ 0.3) between the 
current work and the earlier work. This analysis indicates that with similar experi- 
mental accuracies in the two programs, an uncertainty in temperature between three 
and six times as great should be expected for this program as compared to the previous 
program. Hence, the experimental uncertainties of the two programs are probably of 
the same order . 

An additional uncertainty existed in this program which was not present in the 
earlier program: the sodium line -reversal measurements were made farther from 
the reaction zone than the infrared measurements, far enough so that mixing had begun 
to affect the sodium line temperature. Corrections for this were applied as shown in 
Fig. 5. However, the correction was small (~35°K), so that after the correction was 
applied, the residual error in the sodium reversal measurements was probably signi- 
ficantly smaller than the other uncertainties . 

Thus, the infrared emission -absorption technique for spectral inversion is appli- 
cable to flames as well as to hot gases which exhibit a wide range of temperatures . The 
parametric inversion of infrared spectra should be capable of providing temperature and 
water -vapor profiles as reliable as the general profile shapes that are provided . 

On the other hand, it does not appear that the emission absorption technique in 
the ultraviolet OH band will be quite as useful in flames as the infrared region because, 
within the reaction zone, the electronic states of OH are very substantially out of equi- 
librium with the kinetic temperature of the gas and substantial chemiluminescence is 
occurring [ Refs . 14-18] . Thus, the temperatures measured in the reaction zone by 
emission -absorption in the OH band will be significantly higher than the kinetic temper- 
ature. It is our opinion that, in the current program, a portion of the beam viewed the 
reaction zone, a condition which led to higher emission -absorption temperatures in the 
ultraviolet than in the infrared. If the measurements in the current program had been 
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made somewhat further above the burner face, (where the temperatures would have 
been perhaps 50°K lower) this anomaly would not have appeared. Further experiments 
are required to test this hypothesis . However, the existence of nonequilibrium in the 
OH radicals places fundamental limitations on die regions in which emission -absorption 
probing by use of OH radiation is reliable . In particular, the use of this technique in 
and near combustion zones would not yield reliable results. Unfortunately, the sensing 
of combustion -zone temperatures constitutes an important set of potential uses for such 
techniques, since it is highly desirable to sense temperatures in at least the higher 
temperature regions by short wavelength emission. 

Because of the very rapid variation of the equilibrium OH concentration with 
temperature, a disproportionately large fraction of the OH emission (and absorption) 
originates in the hottest portions of the flame . The result of this condition is that, in 
the inversion of the OH data to obtain temperatures, while the temperature of the hotter 
portions of the flame can be accurately extracted, the extraction of the temperature of 
the cooler portions is subject to very substantial uncertainty. This difference indicates 
that the OH radiation is best suited for the extraction of the maximum temperature of a 
gas having large temperature variations, and the infrared or other techniques should be 
used to extract the temperatures of the cooler regions . 

It is desirable to sense temperatures in at least the higher temperature regions 
of a gas flow from the short wavelength emission primarily because of the greater sen- 
sitivity of the Planck function to temperature at shorter wavelengths . However, if regions 
are to be included in which nonequilibrium of the 1 OH is a possibility, the use of an emission 
absorption technique in the OH band should probably be abandoned. Other techniques in 
the ultraviolet should be investigated as possible means of temperature inversion in non- 
isothermal flames . Temperature inversion based on the relative intensity of the rotational 
lines within a band, has been highly developed for isothermal paths (as used in this report) 
and perhaps could be generalized to yield temperature profiles across nonisothermal paths. 
However, such techniques will sacrifice much of the potential advantage of the. short wave- 
length region, since the absolute emission values, the source of the potentially superior 
accuracy of the ultraviolet region when equilibrium is present, must be discarded where 
nonequilibrium exists. Furthermore, lack of knowledge of accurate transition probabi- 
lities and line widths limit the accuracy of rotational line analyses . It seems that ultra- 
violet line -of -sight probing techniques that are not based upon emission -absorption will 
not be substantially more accurate than infrared measurements and will be at least as 
difficult to perform experimentally. 

Accordingly, future efforts should be directed toward a delineation of those regions 
and conditions in which temperature inversion by OH emission -absorption spectrometry 
is valid (i.e., where OH chemiluminescence is not significant). The various entries of 
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Refs. [14-18] describe observations of various degrees of chemiluminescence, all of 
which indicate OH nonequilibrium, within and near a reaction zone, and then downstream 
decay to complete equilibrium. However, the spatial extent of the chemiluminescence 
varies in the different cases, depending on the reactants and the flame conditions . The 
spatial extent is not known in generaland, as indicated by the current work, could possibly 
be more extensive for diffusion flames than for premixed flames. Thus, it seems that a 
more complete study of the conditions under which substantial chemiluminescence of the 
OH molecule occurs must be carried out. This study could begin by probing an atmos- 
pheric pressure, H 2 _ 02 , diffusion flame of the type used in the present study with both 
sodium and lithium line -reversal measurements and OH emission -absorption measure- 
ments . Such a study would allow the determination of various translational and OH rota- 
tional temperatures and OH ^2 + population throughout the flame. The results would 
indicate the regions in the burner flame where OH chemiluminescence is important and 
hence provide information whereby the currently obtained, high emission -absorption 
temperatures in the ultraviolet could be explained . Further theoretical studies would 
still be required (besides the results in Refs . [14-18] ) to attempt to generalize the con- 
ditions governing OH chemiluminescence. The results of this study would provide a 
means for determining whether a particular physical situation or region would yield OH 
chemiluminescence, and thus a means for deciding whether OH emission -absorption 
measurements with their potentially better temperature accuracy would be worthwhile 
for temperature probing in that region . 

Uncertainty in the temperature determined for a burner segment from the OH 
measurements was strongly dependent on the OH content and hence temperature in that 
segment. The accuracy of the ultraviolet temperature determination is variable along 
the line of sight. An additional analytical study should be performed to determine a 
method by which ultraviolet and infrared emission absorption measurements can be 
combined so that the accuracy of the temperature profile all along the line of sight is 
maximized. 


50 



WILLOW RUN LABORATORIES 


Appendix I 

A BAND MODEL FORMULATION FOR VERY NONUNIFORM PATHS 
1.1. INTRODUCTION 

, The particular band model used for the research described in this report was developed at 
Willow Run Laboratories [1, 2] , and was originally formulated in terms of a nearly-weak or a 
nearly -strong (NW-NS) line approximation. A simple procedure for interpolation between the 
two approximations provided for intermediate optical depths. Applications of this band model 
in its original form yielded anomalous results; at some frequencies, the calculated values of 
apparent radiance with atmospheric absorption exceeded those for the hot source alone. Rea- 
sons for these inconsistencies were determined and corrections made. In the body of the re- 
port, we presented both results of the use of the modified band model and the requirements for 
a more complex interpolation procedure. The purpose of this appendix is to present the details 
of the revised formulation of the model, which represents a closer approximation to the exact 
equation of transfer than the original model, or comparable ones which have been based on the 
Curtis -Godson approximation. 

The rationale in the development of the new model can be illustrated by considering the 
growth of a single, collision -broadened spectral line along a nonisothermal path. Once the 
formulation for such a single line is developed, it is carried over essentially unchanged to the 
case of a random band of overlapping lines. 

The solution of the equation of transfer which describes the radiance of a single, isolated, 
spectral line from a general nonisothermal source under the conditions of local thermodynamic 
equilibrium, negligible scattering, and no significant source of radiation behind the gas, is 

x L 

L^i/) = J X)k (u - v Q , X) exp 

0 

where X is the optical depth (mass per unit cross-sectional area) coordinate along the line 
of sight through the gas. L* (v, X), the blackbody spectral radiance, is a function of frequency, 
or equivalently wave-number, v, and temperature, T = T(X); T in turn, is a function of the posi- 
tion along the line of sight through the gas. The spectral absorption coefficient of the single 
line being considered, k(v - u q , X), is a function of the spectral distance from the line center, 
v - u q , and the optical depth coordinate along the path, X. L ^(v) is the spectral radiance of the 
gaseous source at wave-number u resulting from the single spectral line. X^ is the total opti- 
cal depth of the gaseous source. The total radiance from the line is obtained by integration 
over all frequencies: 


A 

-J k(i/ - v Q , X')dX’ 


dX 


(21) 
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CO X T 


r x 


L = j j* L*(i >, X)k(r - v q , X) exp -j" k(y - u q , X')dX' 

-oo 0 L 0 


dXd(i/ - 


If we introduce the transmittance, defined as 
r X 


t(u - v X) = exp 


-j k(y - V 

_ o 


X') dX' 


(23) 


and interchange the order of integration, assuming that L^(y, X) does not vary significantly 
over the frequency region involved, the equation for L can be written: 


5 


“00 - 

L = - 

CL 

L >o’ X) Jc 

j" t(i> - V Q , X) d (u - U Q ) 

0 

_-0Q ' 


dX 


(24) 


The equivalent width, W, is defined as 

00 

1 


W(X) = 1 - | t(u - u q , X)d(v - v Q ) 


(25) 


-oo 


so that the expression for L becomes 


-/ 

0 


L - 1 L X' x) S dX 


(26) 


Since W is only a function of X, the above equation can be written 


W(X T 


L*(i/ , X)dW(X) 
v o 


(27) 


If we divide the path into N approximately isothermal elements, the expression for L becomes 
N 

Ls VL* [W(X ) -W(X )] .(28) 

/ f V J 1 

i=l 

where X. is the optical path from the observer to the far boundary of the i-th isothermal ele- 
ment and W(X.) is the equivalent width corresponding to that path. 

Calculation of the values of W(X^ shows that, in general, the integrals over X and v - 
are not separable for nonisothermal paths. To allow such separation, an approximation for the 
variation of the line shape along the path is usually made. The Curtis -Godson approximation 
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is most commonly used for this purpose [21, 22, 23] . Another, is the nearly-weak, nearly- 
strong line approximation [1] . The purpose of both approximations is the separation of the 
integrals over- X and v - and the replacement of the nonisothermal path with an equivalent 
isothermal path. Both yield satisfactory values of W(X) for many cases, especially when the 
temperature and concentration variations are not too large. 

In the evaluation of Eq. (28), the accuracy of the integration depends not so much on the 
accuracy of the values of W(X), but on the accuracy of differences between values of W(X), iie., 
on the derivative of W with respect to X. However, an expression which is a reasonable approxi 
mation for W will not always yield a reasonable approximation to dW/dX. To illustrate this, 
we will compute the derivative of W with respect to X, using both the Curtis -Godson approxi- 
mation and the NW-NS approximation. 

First it is necessary to obtain an expression for the equivalent width of a collision -broad- 
ened line in an isothermal path. The absorption coefficient in this case is 


k(u - v Q , X’) 


S(X’) y(X') 

n y 2 (x') + (v - v Q ) 2 


(29) 


where S is the line strength and y its half -width. Upon substitution into Eq. (23) specialized to 
an isothermal path, and with the use of Eq. (25) the following expression for W is obtained after 
some reduction [20] : 


W(X) = 2iry ■ f(x) 


(30) 


where f is the Ladenburg-Reiche function. 

f(x) = xe x [Jq(lx) - iJjUx)] (31) 

SX 

and x = is the dimensionless optical depth. Jq and are Bessel functions of the first kind, 

of order 0 and 1 respectively. Tables of f(x) are found in the literature [25, 26] . 


- 1.2. CURTIS -GODSON APPROXIMATION 

We have used the form of the Curtis -Godson approximation described by Goody [21] . It is 
a method in which the equivalent width of an absorption line as viewed through a nonisothermal 
path can be approximated by the equivalent width through some corresponding isothermal path. 
The approximation involves the definition of an equivalent average half -width and line strength 
for the path extending from the observer to some point, X, in the path. 


(32) 

(equation continued) 


s e (x) 


rx 
J S(X') 
o 


dX' 
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and 


A 

J y(X')S(X) dX' 


r e (» - 


The denominator in the nonisothermal absorption coefficient (Eq. 29), is taken to be independent 
of X', of the form, y 2 (X) + (y - i^) 2 , so that 



X 

J 

0 

1 S(X') dX' 


X 

I 

0 


k(i/ - v , X')dX’ 


y e (X)S e (X)X 


y 2 (x) + {v - y Q ) 2 


(33) 


From this, it follows that 
W(X) = 2;ry e (X)f(x) 

where x, the dimensionless optical depth, is now defined as 
X 

J S(x') dX' 


(34) 


x(X) 


_ 0 _ 


2ny e (X) 


(35) 


Equation (34) can now be used to investigate the behavior of dW/dX in a nonisothermal gas. If 
we differentiate (34) and make use of (32), we get 


dW 

dX 


= S(X) 


df(x) 

dx 


2 - 


y(x) 

r e (» 


f(x)Ty(x) , 

-|y^xr-‘. 


(36) 


The variation with optical depth X can be mapped into a variation with dimensionless optical 
depth, x, and Eq. (36) rewritten as 


1 dW 
S(X) dX 




(37) 


We can now place physical limitations on [1/S(X)] (dW/dX). At X = 0, the quantity [l/S(X)] 
(dW/dX) must be unity corresponding to the region of linear growth for small optical depths. 

It must also be less than one at larger optical depths as the square root region is approached; 
and, finally, it must remain non-negative for all x, so that the equivalent width continues to 
grow as x increases. The variation of [l/S(X)] (dW/dX) in Eq. (37) is shown in Fig. 20 as a 
function of x with y/y g as a parameter. y/y g is the ratio of the half -width at x to its average 
value taken over all the path before it. For y/y g < 2, the function is well behaved, but not for 
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y/y e >2. In the latter case the function becomes larger than one for intermediate values of x, 
indicating a rate of equivalent -width increase greater than that in the linear region.* 


1.2. NEARLY -WEAK, NEARLY-STRONG LINE APPROXIMATION 

This approximation, like the Curtis -Godson approximation, defines an average for the line 
half-width. However, this average is fixed for the path in question and is defined as the average 
over the total path, X : 

1j 


L 

J r (x’)s(x')dx’ 

_ 0 , 

y eL X L 

| S(X')dX’ 

0 

Then the dimensionless optical depth is defined as: 


(38) 


x(X) = 



X 

j S(X')F(X’)dX' 
0 


(39) 


F(X’) is a function which must be unity in the NW approximation and y/y in the NS approxi- 
mation [5] . A form must be chosen for the variation of F(X') so that it produces a reasonable 
transition between the two approximations for intermediate optical depths. The equivalent 
width is given again as 

W(X) = 27ry eL fW (40) 


If we take the derivative of Eq. (40 for the NW approximation we obtain 

g = S(X)« (41) 

dX dx 

and for the NS approximation, we obtain 

dw s(x) r(x)dM (42) 

dx r eL dx 


The values for [l/S(X)] (dW/dX) (according to Eqs. 41 and 42) are shown in Fig. 21. It can be 
seen that the NW approximation never becomes physically unrealistic, i.e., [ 1/ S(X)] (dW/dX) > 1 


” *These errors, mentioned by Goody [2T] , are treated in detail by Drayson [22] . However 
their analyses dealt only with the absorption (or equivalent width). On the other hand. Walshaw 
and Rodgers [26] analyzed the effect of the Curtis -Godson approximation on the derivative of 
transmittance with optical depth, for several band models. 
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DIMENSIONLESS OPTICAL DEPTH, x(X) 

FIGURE 20. EQUIVALENT -WIDTH DERIVATIVE IN THE 
CURTIS -GODSON APPROXIMATION 
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FIGURE 21. EQUIVALENT -WIDTH DERIVATIVE IN THE NEARLY -WEAK 
AND NEARLY-STRONG LINE APPROXIMATIONS 
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however,' for the NS approximation, the quantity [l/S(X)] (dW/dX) becomes greater than one at 

small values of x for — — > 1. This is accounted for by the proper choice of F(X'). However, 

, y eL 

we are left with no information upon which to base F(X') at this stage. 

It appears that neither the NW-NS approximation (unless F(X) is accurately specified) nor 
the Curtis -Godson approximation give sufficiently accurate values for dW/dX over the full range 
of optical depths and for all possible values of y/y . On the other hand, it should be noted 
that the anomalies and. errors in the curves for (l/S)(dW/dX) as given by the Curtis -Godson 
approximation are not serious as long as y/y never exceeds a value of about two. This con- 
dition exists [ 21 , 22 , 27 ] for most atmospheric paths but not for highly nonhomogeneous paths. 
For example, hot products of combustion viewed through long atmospheric paths yield values 
of y/y g as large as 20 (Appendix III). 

To obtain a more acceptable expression for dW/dX, consider the growth of an isolated 
spectral line. The equivalent width is given exactly by 



CO 

r x -i 

r 

W(X) = 1 - 

exp 

- k(Ay, X’) dX' 
J 

J 

30 

L o J 


d(Ay) 

and its derivative, in terms of a Lorentz profile, by 


(43) 


00 

r 


dW 

dX 


S(X)y(X) 
y 2 (X) + (Ay 2 ) 


exp 


-00 


-A 

M 

L 0 


S(X')y(X') 
y 2 (X') + (Ay) 2 


dX' 


d (Ay) 


(44) 


where Ay = y - y . In the original treatment of the Curtis -Godson and the NW-NS approxima- 
tions, the separation of frequency -dependent and path -dependent variables was accomplished by 
substitution of an effective value, y , in both parts of the integrand. The only difference in our 
treatment from previous treatments is that we will use such an approximation only in the ex- 
ponential term. Thus, in terms of y g = y g (X) defined by Eq. (32) and x defined by Eq. (35), 

Eq. (44) becomes 

CO 


dW 

dX 


if 


exp 

y \X) + (a vY 


y (X) 

-2 xyfflj-S 2 

yf(X) + (Ay^ 


d (Ay) 


(45) 


-CO 

This is equivalent to making a Curtis -Godson substitution in the expression for the derivative 
of W, rather than in the expression for W itself. The process of taking a derivative after an 
approximation has been made is reversed, a treatment which should be inherently better than 
the techniques described earlier for the calculation of radiances. This new approximation can 
be compared to the previous two if dW/dX is considered as the limiting value of the contribution 
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to W from a small element of path, AX, viewed through the path ahead of it. In the previous 
two approximations, the line shape of the element AX and the path in front of it are both given 
by the line shape for the equivalent, homogeneous path. In the current approximation. we as- 
sumed that the element AX has its true line shape, while the path in front of it has an equivalent 
homogeneous shape. Equation (45), checked against the exact relation (Eq. 44), for a two-layer, 
nonisothermal path, has been found to be very accurate. Further investigation of its accuracy 
is to be the subject of future work, It would be possible to fully develop a procedure for treat- 
ing isolated Lorentz lines. However, our interest is ultimately in a random band model which 
treats the average of many lines. 

1.3. BAND-MODEL CONSIDERATIONS 

The calculation of an average spectral radiance at the boundary of a radiating body of gas, 
in a frequency interval encompassing a number of rotational lines, is represented by the rela- 
tion 


co 

L, = Jgfc’ - u)h u (u')dv' (46) 

0 

where g (u, u') is an averaging function (e.g., the slit function of a spectrometer of moderate 
resolution, such as those which produce original absorption spectra from which band-model 
parameters are extracted). The exact spectral radiance, L^, is given by the equation of trans- 
fer, Eqs. (21) and (23). Thus, Eq. (46) can be written 

« r(X L ) 

L^ = J g(u' - v) J L*(t/, r’JdT'di/ (47) 

0 1 

where r is the perfectly resolved transmittance, and the primes denote variables of integration. 

The order of integration in Eq. (47) can be inverted, and when frequency intervals are small 
enough that the Planck function is essentially constant, the result is the band-model expression: 

f « L > 

L (v) = L* (u, t') dr ' (48) 

v J u 

1 

where 

00 

t(v) r J g(v' - v)t{u')Au' (49) 

0 
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Thus, r is seen to be the spectral transmittance measured by a conventional instrument or cal- 
culated by a band -model. 

Equation (49) can be evaluated to obtain the average spectral radiance, L^, by first comput- 
ing r at a number of points along the path. The blackbody radiance is then represented as a 
function of r, and integrated numerically over t, with, for instance, the use of a simple sum- 
mation. The various band models in current use, including those based on the Curtis -Godson 
approximation and the present model in its original form, do yield values for the transmittance 
with an accuracy sufficient for most engineering applications. However, as for the case of isolated 
lines, Eq. (49) involves a differentiation of the expression for r with respect to the optical path. 
Differentiation, of course, accentuates errors and uncertainties; in the isolated line case pre- 
sented above, the derivative of the equivalent width exhibited the error. Therefore, an expression 
which yields reasonable values for r will not necessarily yield reasonable values for dr and, 
hence, for L^. It was necessary, therefore, to develop an expression for r which has a derivative 
which will yield more realistic values of the difference in transmittance with a specified in- 
cremental increase in physical path. 

The transmittance for a band of spectral lines with some overlapping can be expressed as 

[2ll 

r = exp (r n - 1) (50) 


where r^ is the transmittance for an equivalent band of well-isolated nonoverlapping lines. 
Differentiation of Eq. (50) with respect to the optical depth, X, yields 


dr 

dr _ - n 
dX ~ T dX 


(51) 


which indicates that, with a reasonably accurate value for r, the accuracy of dr/dX depends 
directly on the accuracy of dr^/ dX. 

In our current work, we have used the NW-NS approximation exclusively and at this point 
our consideration is limited to this model to find specification of the function F(X') which yields 
accurate values of a transmittance derivative. For the NW-NS approximation [5] 

f n = ^eL f( ^ (52) 

where the argument of the Ladenburg-Reiche function :is,. in this case , 


x(X) 


tl 


k(X') 


g(X') 

L S eL . 


]7?(X’) 


dX' 


(53) 
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in which k = k(i/, T) is the first band-model parameter, the average absorption coefficient, 
identified as the average line strength to spacing ratio, and 0 = p(v, T) is the second band- 
model parameter, the line-overlap factor, identified as 2tt times the average line width to spac- 
ing ratio. The effective value 0 e ^ i s defined analogously to y for the case of a single line: 


0 


eL 


j k(X')/3(X') dX' 
0 


J k(X') dX' 
0 


(54) 


77 (XM 

F(X') has here been replaced by [0(X')/0 eL ] ” , a suitable functional form for F(X'). Here 77 

is an interpolation parameter between the nearly weak and nearly strong approximations, varying 
from 0 for the nearly weak to 1 for the nearly strong. 

In an earlier study [5], an investigation of the form for 77 is described. As a result of that 
study, several reasonably useful empirical forms for 77 were determined. However, the infor- 
mation used to make these determinations was not sensitive enough to values of 77 to provide 
accurate information about its form. One of these empirical forms is: 


77 = 0 for x < 1 

77 = (x - l)/9 for 1 s x s 10 ' 

77 = 1 for x > 10 

This appears to make 77 a function of x rather than X', and the determination of x appears as 
an iterative procedure. However, in practice the integral for x is evaluated by a summation, 
and a value for 77 for a particular term is determined from the value for x calculated from the 
sum of the preceding terms. Thus, 77 is uniquely defined for each X'. 

Differentiation of Eq. (52) and utilization of Eq. (54) yields 


dr 

n 

dX 


df_ 

dx 


'Mi 


eL 


M) 


k(x) 


(55) 


which shows that the derivative of dr^/dX varies with the absorption coefficient (as would be 
expected), the derivative of the function f(x), and the local values of 0(X) and tj(X). The phys- 
ical significance of Eq. (55) can be seen in Fig. 22, a normalized plot of -(l/k)dr /dX versus 
x, with 0/0 eL as a parameter. In the limit, as x — 0, the band model reduces to a linear growth 
law, so that -(1/k) dj^/dX approaches unity. This function must be less than unity for all values 
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DIMENSIONLESS OPTICAL DEPTH, x(X) 

FIGURE 22. TRANSMITTANCE DERIVATIVE IN THE 
ORIGINAL BAND MODEL 
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of X greater than zero, as the square root growth law becomes operative. However, Fig. 3 shows 
that the empirical specifications for r) do not conform to this requirement for certain values of 
an d x > an inconsistency analogous to that shown earlier for isolated line growth, dr^/dx 
is correctly given by [21] 


dr 

jn 

dX 


1 dW 
5 ' dX 


(56) 


where W is an average equivalent width and 6 is an average line spacing in the spectral interval 
under consideration. For lines of equal strength and width, after some rearrangement, Eqs. 

(56) and (45) yield 


dr 
n 

dX 


00 

r 


k(X)j3(X) 

o o exp 

8 (X) + (2rAiV6r 


-00 


-2x0 e (X) 


(X) 


1 


£ 2 (X) + (2ttW<5) 2 


(2ttAu 
d \— 


(57) 


where x denotes the dimensionless optical depth as defined in Eq. (35) and used in Eq. (45). 
Rewritten in terms of the two bandmodel parameters, 


X 

| k(X') dX’ 


_ 0 

= ^e (X) 


(58) 


where /3 g (X) is the average value of the second band-model parameter over the path from the 
observer to the point X, similar to the definition of 7 e (X) in Eq. (32) 

X 

jk(X')|3(X')dX' 

/3 e (X) = 9-^ (59) ' 

J k(X') dX' 

0 


The distinction between (3^ (a constant determined over the whole path), defined in Eq. (54) and 
(3 e (X) (a function of X, the optical depth coordinate) should be noted. 


The change of variables tan 8/2 = 2n(Ai//p e 6) and a trigonometric substitution yields 


1 dT n _ 2 J_ 

r 

exp f-x(l + cos 9] 


k(X) ' dX ^ P e 

J 

f 

1 \f) 2 ^} 

L^ e 

+ 

m - '} 

cos 9 


-d0 


(60) 
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Equation (60) has been integrated numerically; the results are plotted in Fig. 23, as. -(l/lOdr^/dX 
versus x; 0 (X)/B^(X) is the parameter. The resulting curves exhibit none of the anomalies 
seen in Fig. 22; the curve for = 1 is identical in both figures. 

The interpolation parameter, rj, can now be evaluated by equating Eqs. (60) and (55). These 
two equations involve different sets of variables, i.e., Eq. (60) is given in terms of x and 0 e (X), 
and Eq. (55) is given in terms of x and <3^. The relationship between x and x is complex and 
depends on ~q, /3 g (X) and /3 gL in an integral form. Furthermore, although x is a unique function 
of the optical path extending from 0 to X, x is not a unique function in the same sense, since it 
depends on the path between X and X^ through /3 gL as well. Such dependence is not physically 
realistic but appears as an integral part of the NW-NS approximation. In determining a form for 
f], this dependence must be accounted for so that the computed value of (l/k)(dT n /dX) is unique 
despite the nonunique character of x, although this is very difficult to insure, except for certain 
types of paths. In particular, if 0 differs significantly from its average value, 0 . , in only a 
very small portion of the optical path, then the contribution from that portion of the path to the 
integral defining x is negligible compared to the contributions from the remainder of the path. 

In this case, the following approximate relation holds: 


X 


k(X’) 


■ §0c)f (x ' dX 

[_0eL J dX 


A 

0 


k(X') dX' 


The above restriction also implies that 

0 (X) = constant = 0 T 
e 'eL 

and that 

x - x 


Thus, if 0 (X) is not a strong function of X and if large variations in 0(X) are allowed in very 
small fractions of the optical path only, values for tj can be found in terms of x and 0/0 e ^ (ex- 
cept for 0/0 = 1 where ?? is indeterminate). The resultant values of r\ are shown in Fig. 24. 
Use of these values in Eq. (55) will insure that the derivative dr/dX will equal that given by 
Eq. (60). For convenience in application, we made simple analytical approximations to the 
curves of Fig. 24 by fitting them with functions as follows: 


x + 

lj(x, 0/0 ) = — 

e X + 


Vgggej 

A 2 (fi/|3 e ) 


(61) 
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FIGURE 23. TRANSMITTANCE DERIVATIVE IN THE 
MODIFIED BAND MODEL 



DIMENSIONLESS OPTICAL DEPTH, x(X) s x (X) 


FIGURE 24. INTERPOLATION PARAMETER AS A FUNCTION 
OF OPTICAL DEPTH AND LINE-OVERLAP FACTOR 
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^A 3 (0/0 e ) 

r)(x, 0/0 e ) = A 5 (/3//3 e ) + xA 4 (0/0. e ) 

For any particular instance, the choice between Eqs. (61) and (62) is made on the basis of the 
quantity 

A g (|B/0 e ) = 1 +0.185 (3/p e (63) 

If x > E, Eq. (61) is used; otherwise Eq. (62) is used. Values for A^, A 2 , Ag, and A' 4 are tab- 
ulated in Table 11. Values of dr n /dX obtained when the values in Table 11 are used in place of 
the values shown for t) in Fig. 5 are sufficiently accurate for most applications of this band- 
model. 

The limitation on the variation of 0(X) and 0 e (X) mentioned earlier is not as severe as it 
might first appear. For instance, in cases where a relatively small hot gaseous source is to 
be viewed through a long, cool path, the value of (3 is almost invariant over the whole atmo- 
spheric path which contains the largest part of the optical path. Only in the relatively small 
portion of the path comprising the hot gaseous source does /? vary substantially. 

The introduction of an interpolation parameter has little physical basis and is not the only 
way Eq. (60) could be introduced into a computational procedure. An alternative method for use 
of the above results would be the direct application of Eq. (60) into a band-model computer code. 
The resultant values of dr^/dX could then be used in Eq. (51), if T n in Eq. (50) is evaluated in 
terms of x and /3 g (X) . Such a procedure would eliminate the need for a distinction between the 
NW and NS approximations and, therefore, eliminate an interpolation procedure. Unfortunately, 
when Eq. (60) is used, an analytic form for the solution of the integral is not easily obtained. 
Although a solution has been obtained in the form of an infinite series, this represents little 
improvement over a table of numerical values. Whether there would be a net gain in computa- 
tional efficiency remains to be seen. However, a direct application of Eq. (60) might be ex- 
pected to be somewhat more satisfactory because the' restriction on the variation of (3 g (X) com- 
pared to /3 gL required in the evaluation of rj would not be needed. This alternative mode of the 
band-model formulation is one of the subjects of a continuing investigation. 

In summary, we have shown how both the Curtis -Godson approximation and the NW-NS 
approximation can introduce anomalies into the radiance computed over a highly nonisothermal 
path for an isolated spectral line. The application of a Curtis -Godson approximation to the ex- 
pression for the equivalent -width gradient along the path produced a formulation yielding greatly 
improved radiance values which exhibit no anomalies. In the generalization to a random band 
model, the expression concerning the equivalent -width gradient has been used to obtain an inter- 
polation between the nearly -weak and nearly -strong approximations valid for that band model. 
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The most important result presented herein is the form of the transmittance gradient 
given in Eq. (60). It is the basis of a computational procedure for computing the radiance from 
highly nonisothermal gaseous sources. This procedure yields significant improvements over 
previously used techniques. It has been shown to be a better approximation than those more 
commonly used and yields physically realistic results for all values of the path and band-model 
parameters. 


TABLE 11. COEFFICIENTS TO BE USED IN 
EQS. (61) AND (62) 


0/0 e 

A 1 

A 2 

A 3 

A 4 

0 

-0.68 

-0.35 

0.55 

0.07 

2 

-0.37 

0.62 

0.71 

0.43 

5 

0.60 

3.15 

0.90 

0.79 

10 

3.19 

9.25 

1.15 

1.30 

20 

9.10 

22.8 

1.63 

2.26 

50 

29.1 

67.1 

2.78 

4.56 

100 

105 

209 

4.05 

7.09 
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Appendix II 

COMPARISON OF THE MODIFIED BAND MODEL WITH OBSERVED HOT GAS SPECTRA 

The measured spectra in Figs. 25 through 30 are those of furnace -heated gas samples, 
60-cm long, obtained in the previous study [1] . We obtained the band-model predictions (the 
less detailed spectra in each figure) using Eqs. (1) and (2), modified as indicated in Appendix I. 
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FIGURE 25. SPECTRAL RADIANCE, RUN 10216807. 
Gas: HgO/COg/Ng; P = 57, 28, 675 mm; temperature 
(OK) = 378, 537, 723, 958, 1127, 1158, 990, 752, 555, 383. 
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FREQUENCY (cm" 1 ) 

FIGURE 26. SPECTRAL ABSORPTANCE, RUN 10216807. ' 
Gas: H 2 0/C0 2 /N 2 ; P = 57, 28, 675 mm; temperature 
(°K) = 378, 537, 723, 958, 1127, 1158, 990, 752, 555, 383. 



FIGURE 27. SPECTRAL RADIANCE, RUN 10216811. 
Gas: HgO; P = 670 mm; temperature (°K) = 382, 537, 
723, 953, 1128, 1160, 990, 751, 558, 389. 
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FIGURE 29. SPECTRAL RADIANCE, RUN 10216812. FIGURE 30. SPECTRAL ABSORPTANCE, RUN 10216812. 

Gas: COgl P = 760 mm; temperature (°K) = 386, 528, Gas: COg: P =760 mm; temperature (°K) = 386, 528, 719, 

719, 953, 1130, 1160, 979, 737, 541, 387. 953, 1130, 1160, 979, 737, 541, 387. 
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Appendix III 

CALCULATION OF RADIATION FROM HOT H 2 0 AND C0 2 VIEWED THROUGH 
A COOL INTERVENING ATMOSPHERE 

There are a number of engineering problems requiring quantitative calculations of the 
thermal emission from hot water vapor and carbon dioxide, viewed through atmospheric paths 
in which these species also appear as absorbers. In recent years, molecular band models, 
originally developed for atmospheric transmission calculations, have been extended to handle 
the case of inhomogeneous bodies of hot gases which emit as well as absorb. This appendix is 
concerned with their further extension to the case of a hot body of gas imbedded in a cool atmo- 
sphere, treated as a single radiative -transfer problem. The computational method is illustrated 
by a comparison of predicted and. observed spectra in the 2. 7 -pm region for a small methane- 
oxygen flame. • 

The arrangements for the illustrative experiment were as follows: the flames were gen- 
erated with a multiple-diffusion burner originally built for hydrogen -fluorine combustion [11]; 
it operated equally well with methane and oxygen. The burner was installed in the sample beam 
of the double-beam spectrometer described in [1]. An auxiliary system of mirrors, two diagonals 
and a sphere, provided for measurements of spectral radiance and absorptance of the flame 
viewed directly and through a 10.32-m path of ambient air; with the flame extinguished, the ab- 
sorptance of the atmospheric path alone could also be measured. The observed spectral radiance 
in the 2.7-pm region of a stoichiometric CH^/0 2 flame with and without an intervening atmo- 
spheric absorption are shown in Fig. 31. The extreme wings of the flame emission band, com- 
prised primarily of lines originating from higher rotational levels, are relatively unattenuated, 
in contrast with the central region in which the low temperature atmospheric absorption is en- 
countered. This characteristic shape for a a vibration-rotation band of a highly nonisothermal 
path is perfectly analogous to the familiar shape of a self -reversed spectral line under corre- 
sponding conditions. 

For the band-model calculations, a flame model, similar to that used in [11] , consisted of 
an isothermal core and a mixing region in which temperature and the concentration of the com- 
bustion products decreased parabolically, as the result of dilution with ambient air, as shown in 
Fig. 32. However, in this case, recombination of two dissociated species, HgO and COg, had to 
be accounted for. In order to avoid a detailed equilibrium calculation for all species concen- 
trations at each local value of temperature in the mixing region, the assumption was made that 
the reactions occur independently: 


1/2H 2 +OH; H 2 0 

(64) 

i/20 2 + co r co 2 

(65) 
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FIGURE 31. OBSERVED AND PREDICTED SPECTRA FOR A 2.5-cm 
DIAMETER CH 4 /0 2 FLAME VIEWED DIRECTLY AND THROUGH A 
10-m PATH. The smoother curves are the band-model predictions; 
the upper curve is the flame alone; the lower curve is the flame and 
atmospheric path; the middle is the product of flame radiance and at- 
mospheric transmittance which have been calculated separately. 



FIGURE 32. MODEL FOR FLAME STRUCTURE 
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This permitted the approximations P(OH) ~ 2P(H2) and P(CO) = 2P(0 2 ) which, combined with 
the definitions for the equilibrium constants and the over -all reaction equation, provided the 
relations 

P(H 2 0) = 0.667 - 3[K 1 P(H 2 0)/2] 2//3 (66) 

P(C0 2 ) = 0.333 - 3[K 2 P(C0 2 )/2] 2/3 (67) 

Equations (66) and (67) were solved numerically at each temperature, with values for the equi- 
librium constants, K = K(T), given in [27] . 

As inputs to the flame model, the core temperature, T , was taken as that measured in the 
center of a comparable CH 4 /0 2 flame by a sodium line-reversal technique; the core radius 
r c /R was evaluated in accordance with the well-known properties of turbulent mixing of jets 
discharging into a quiescent medium [28] • The resultant values, T^ = 2800°K and r^/R = 0.74, 
were then used in the band-model calculations of the spectral radiance of the flame viewed 
directly and through the atmospheric path. The resultant band-model predictions, shown in 
Fig. 31 as the upper and lower smooth curves, compare favorably with the observed spectra. 

Also shown in this figure, the middle smooth curve is the product of the flame radiance and 
atmospheric transmittance calculated separately. This procedure is incorrect in principle, since 
an atmospheric transmittance factor, calculated with such a band-model formulation (or for that 
matter determined experimentally by conventional techniques), refers to the attenuation of a 
spectrally continuous source. In this case, the hot gaseous source is spectrally discontinuous, 
with a high degree of correlation in spectral line content with that of the absorbing atmosphere. 
Hence, except for spectral regions in which lines are so close as to completely overlap, the 
application of an atmospheric transmittance factor, so determined, results in an overestimate 
of the radiance at a distant point. For the relatively short atmospheric path in the example 
shown, the two methods yield results that do agree to within the experimental error of the mea- 
surements. However, these effects become more pronounced at longer ranges. To illustrate, 
a set of calculations was made for the spectral radiance of a hypothetical source 1 m in diam- 
eter with a parabolic profile varying from 1200°K in the center to 300°K at the boundary, and a 
corresponding composition profile based on a reasonable mixture of combustion products of 
methane in air, observed at ranges of 10, 100, and 1000 m through a standard sea-level atmo- 
sphere at 50% relative humidity. The results are shown in Fig. 33, in which the upper curves 
are the radiance spectra of the hot source with no atmospheric absorption; the lower and 
middle curves are the apparent radiances at the indicated ranges, respectively calculated 
by the correct and incorrect procedure. It is clear that in selected spectral bandpasses, 
the apparent radiance would be lower by an order of magnitude or more than the product of 
the actual radiance and atmospheric transmittance. 
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FIGURE 33. APPARENT SPECTRAL RADIANCE OF A HOT SOURCE 
VIEWED THROUGH VARIOUS ATMOSPHERIC PATHS. The upper 
curves represent the hot source alone; the lower curves represent the 
combined hot source and atmospheric path; the middle curves show the 
product of the hot source radiance and atmospheric transmittance which 
have been calculated separately. 
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This appendix indicates the importance of the proper treatment of the apparent radiance 
distribution or spectral radiant intensities of distant hot combustion -product sources viewed 
through the atmosphere. The hot -gas band model developed under this and preceding contracts 
does provide a basis for such treatment, and is currently in use in various studies of rocket and 
jet-engine exhaust radiation. 
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Appendix IV 

BAND -MODEL PARAMETERS FOR LOW TEMPERATURES 

The values of the band-model parameters in this study were extracted from the General 
Dynamics tabulations [4]. These parameters were determined from the analysis of spectral 
absorptances of isothermal samples of hot gases and are reasonably reliable throughout the 
temperature range of the original measurements. However, the tabulations were extended down 
to 300°K by extrapolation. When these values were used in the present band model for the esti- 
mation of transmittances in atmospheric paths at some wavelengths, the results were in con- 
siderable disagreement with calculations made with the atmospheric transmission models and 
parameters described by Anding [29] • 

Since the latter have been demonstrated to be reliable over large ranges of atmospheric 
paths, horizontal and slant, the parameters in the hot gas band model were adjusted until a sat- 
isfactory agreement was achieved. Figure 34 shows the results of this adjustment for a 29.2 -km 
slant path through a standard atmosphere; the path parameters are listed in Table 12. Figure 35 
shows a comparison of the two models applied to a path which has the same temperature and 
pressure profiles, but which is shorter by a factor of 100. 


TABLE 12. PATH PARAMETERS IN THE HOT GAS 
BAND MODEL 


Distance 

Temperature 

p(h 2 o) 

p(co 2 ) 

p(n 2 ) 

(cm) 

(°K) 

(atm) 

(atm) 

(atm) 

0 

288.0 

0.956E-2 

0.32E-3 

1.0 

2.92E5 

281.7 

0.676E-2 

0.277E-3 

0.886 

5.86E5 

275.2 

0.399E-2 

0.252E-3 

0.785 

8.7E5 

268.7 

0.244E-2 

0.222E-3 

0.693 

1.17E6 

266.2* ' 

0.155E-2 

0.195E-3 

0.609 

1.462E6 . 

255.7 

0.612E-3 

0.171E-3 

0.535 

2.93E6 

223.3 

0.117E-4 

0.063E-3 

0.199 
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FIGURE 34. TRANSMITTANCES FOR A 29. 3 -km ATMOSPHERIC PATH 



1800 2000 2200 2400 2600 2800 3000 3200 3400 3600 3800 4000 


v (cm *) 

FIGURE 35. TRANSMITTANCES FOR A 0.293-km ATMOSPHERIC PATH 
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Appendix V 

ULTRAVIOLET EMISSION AND ABSORPTION SPECTRA 

In Figs. 36 through 40 ultraviolet emission and absorption spectra are shown that were obtained 
for the 5 burner profiles (shown in Figs. 10 through 15). On all absorption curves, the absorp- 
tion zero is determined from the local position of the trace between lines, and full absorption 
is taken to be the lower limit of the chart. On emission spectra, emission zero is the lower 
limit of the chart. The radiance calibration curve shown on each emission trace is the trace ob- 
tained from the corresponding observation of the spectral radiance standard lamp. The calibra- 
tion curve shown must be divided by the calibration gain factor shown on each emission spectrum 
in order for the emission spectrum and the calibration curves to be comparable. (The radiance 
of the standard lamp was always much less than the radiance of the flame in the OH lines.) The 
radiance values of the standard lamp as determined by NBS, interpolated to the frequencies of the 
seven OH lines considered are given in Table 13. In the data reduction, these values were re- 
duced by a factor of 0.92 to account for reflection losses in the lens between the lamp and the 
measurement system. 


TABLE 13. SPECTRAL RADIANCE VALUES FOR THE 
SPECTRAL RADIANCE STANDARD LAMP 


Line 

Frequency 

(cm -1 ) 

Spectral Radiance 
(pW/ cm^-sr-cm' 1 ) 

V 

32415.5 

0.349 

R 2 3 

32489.5 

0.337 

R 2 4 

32517.6 

0.333 

R 2 7 

32572.6 

0.324 

R 2 14 

32551.7 

0.327 

R 2 15 

32531.8 

0.331 

R 2 20 

32363.5 

<*•? 

0.358 
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